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*0 ' A framework allowing for perturbative calculations to be carried out for quantum field 

theories with arbitrary smoothly curved boundaries is described. It is based on an expan- 
0\ ', sion of the Green function for second order differential operators valid in the neighbourhood 

_, , of the boundary and which is obtained from a corresponding expansion of of the associated 

r^ ' heat kernel derived earlier for arbitrary mixed Dirichlet and Neumann boundary condi- 

tions. The first few leading terms in the expansion are sufficient to calculate all additional 
divergences present in a perturbative loop expansion as a consequence of the presence of 



O . the boundary. The method is applied to a general renormalisable scalar field theory in four 

^ ' dimensions using dimensional regularisation to two loops and expanding about arbitrary 

background fields. Detailed results are also specialised to an 0{n) symmetric model with 
^ ' a single coupling constant. Extra boundary terms are introduced into the action which 

give rise to either Dirichlet or generalised Neumann boundary conditions for the quantum 
fields. For plane boundaries the resulting renormalisation group functions are in accord 
with earlier results but here the additional terms depending on the extrinsic curvature of 
the boundary are found. Various consistency relations are also checked and the implica- 
tions of conformal invariance at the critical point where the (3 function vanishes are also 
derived. For a general scalar field theory, where the fields (p attain specified values ip on the 
boundary, the local Scrodinger equation for the wave functional defined by the functional 
integral under deformations of the boundary is also verified to two loops. The perturbative 
expansion for the wave functional is defined by expansion around the solution of the classi- 
cal field equations satisfying the required boundary values and the counterterms necessary 
to derive a finite hamiltonian operator, which includes a functional laplace operator on the 
fields (^, are found to the order considered. Consistency of the local Schrodinger equation 
with the renormalisation group to all orders in perturbation theory is also discussed. 



1 Introduction 

Although quantum field theory has an enormous literature, papers devoted to the 
properties of and perturbative calculations for quantum fields in the vicinity of a boundary 
are relatively sparse. Nevertheless considering quantum field theory on a manifold with 
boundary is not without interest for the following reasons: 

(i) The Casimir energy [1,2,3] is obviously an effect which depends on the presence of 
boundaries and exemplifies the crucial zero point energy of quantum fields. However 
explicit calculation for non trivial geometries and for interactions are difficult and 
the usual renormalisation framework requires extension. 

(ii) In statistical physics there are additional critical indices associated with phase 
transitions for physical observables measured near a boundary [4]. These may be 
calculated in the framework of the e expansion of the same quantum field theory 
used to obtain the bulk exponents but in the presence of a boundary [5] . 

(iii) In bag models, which may possibly be regarded as an approximation to QCD with 
some phenomenological validity, it is necessary to consider quantum fields inside 
cavities with suitable boundary conditions [6] . Various perturbative calculations in 
bag models have been undertaken based on using a multiple scattering expansion 
for the propagator in the presence of a boundary [7] . 

(iv) Open strings correspond to conformal field theories with additional operators at- 
tached to the boundary of the two dimensional world sheet [8] . 

(v) Cardy has shown how additional information on the structure of two dimensional 
conformal field theory may be obtained by considering the usual machinery of 
operator product expansions with boundary operators [9]. 

(vi) An alternative formulation of quantum field theory to the conventional approach 
may be given in terms of a Schrodinger like representation with wave function- 
als \E'(<^) defined by the functional integral over quantum fields on a manifold M. 
which attain a specified value (p on the boundary dM.. Such wave functionals obey 
formally a Schrodinger like equation where the time differentiation corresponds 
to deformations of the boundary [10,11]. However the second order functional 
derivatives in the Hamiltonian operator involve singularities which require careful 
discussion. The solutions of the functional Schrodinger equation provided by the 
functional integral over fields on manifolds with a boundary are essentially identical 
to formal constructions of wave functionals satisfying the Wheeler-DeWitt equation 
in quantum gravity [12] which have been extensively explored recently. 

However, in the presence of a boundary the usual efficient calculational techniques 
based on momentum space representations of the propagators are no longer valid, although 
in the case of plane boundaries there are possible extensions [5] . Here we wish to discuss 
quantum field theory on curved space with an arbitrary smoothly curved boundary. For 
simplicity at this stage we restrict our attention to a general renormalisable scalar field 
theory on a Euclidean space with a positive definite metric in four dimensions using di- 
mensional regularisation, developing further the treatment some time ago of Symanzik 



[10]. However, our methods should be feasible for more interesting field theories, such as 
a models in two dimensions which are relevant for strings. 

In order to carry out calculations for the short distance behaviour of amplitudes for 
a quantum field theory in the neighbourhood of a boundary it is crucial to be able to 
determine besides the structure of the leading singularities at coincident points x' ^ x 
of the Green function Ga{x,x') for the elliptic operator A, defined when the action is 
expanded to quadratic order around some background, also the behaviour as x, x' approach 
the boundary. Recently we have developed [13] an alternative to the multiple scattering 
expansion based on an analysis of the heat kernel GAiXjX';T) corresponding to e"^"^ as 
r ^ in the neighbourhood of a boundary for either Dirichlet and Neumann boundary 
conditions. The results depend on an extension of the usual DeWitt ansatz [14] to include 
an effectively semi-classical expansion about geodesies which undergo reflection at the 
boundary. The details of this expansion are summarised in appendix A. Throughout we 
maintain manifest reparameterisation invariance with respect to the coordinates x^ for dAi 
and for the purposes of calculation, in the neighbourhood of dA4, adopt the coordinates 
x^ = {x\ y) where y is the geodesic distance from x^ E A4 to x^ = x^ E dM. defined by 
the geodesic tangent to the unit normal n^ at a;\ In these coordinates the metric g^i, on 
M. is given by 

ds"^ = 7,^(x, y)dx'dx^ + dy"^ , n^ = (0, 1) , (1.1) 

where 7ij(x, 0) = 7ij(x) is thus the induced metric on dM. 

We assume that the operator A, acting on sections of some vector bundle over M., 
is of the form 

A = -1)2 ^ X , D^ = -^D^g^'-^D^ , D^ = d^ + A^, (1.2) 

yd 

for X, A^ matrix valued fields on M.. A^ can be regarded as an external background 
gauge field but may also be taken to include any appropriate spin connection. Assuming 
also that A is a symmetric operator acting on vector fields ^{x) over M. the general form 
of the boundary conditions considered here, which is sufficient for most field theoretic 
applications, involve ^ and its normal derivative on dM. in the form 

(Pn^I)^+V)e|9_^=0, (1-P)^|^_^=0, V^ = V, Pij = i;P = i;, (1.3) 

for 'P(x) a projection operator and 'i/'(^) a matrix valued function on dAi. Clearly from 
(1.3) taking V = corresponds to pure Dirichlet boundary conditions while P = 1 is the 
Neumann case. With these boundary conditions the Green function Ga is defined by 

A^Ga(x, x') = 5^(x, x') = —6'^(x - x') . (1.4) 

In the interior of Ai the DeWitt asymptotic expansion of the heat kernel Qa deter- 
mines a corresponding expansion of Ga which for discussing renormalisable field theories 
in four dimensions may be restricted to the form [15] 

2 

GA = J2Gna^ + GA , (1.5) 

n=0 



for a^{x,x') the Seeley-DeWitt coefficients, nonsingular for x ~ x', and the singular part 
of Ga is given entirely in terms of Gn{x,x') where 
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for a"(a;,x') the geodetic interval given by i J^ ds x{s)'^ for a;'^(s) the geodesic path from 

x^(0) = x' to x^(l) = X. A2(x,x') is a nonsingular symmetric biscalar which may be 
absorbed in a^ but the expressions (1.6) are convenient for correspondence with previous 
results. The representation provided by (1.5,6) is valid for an arbitrary dimension d oi Ai, 
with dAi of dimension d — 1. However in the expressions for Gi, G2 the pole in the first 
term of each on the r.h.s. of (1.6) at e = 4 — d ^ has been subtracted so that there 
is a well defined limit for d = 4, this ensures a convenient form for Ga which is regular 
for £ — !► as well as a; ~ x'. The subtraction also represents the removal of infra-red 
divergences of the corresponding flat space Fourier transforms /c~^, k~^ when k ^ 0. In 
these pole terms an arbitrary mass scale fu. has been introduced for dimensional consistency, 
as usual in dimensional regularisation. Ga is independent of /U, any variation in fj, may be 
compensated by a corresponding change in Ga in (1.5). 

In discussing a renormalisable field theory it is also sufficient to resort to just the 
first few terms in a covariant Taylor expansion of an{x,x') about x' = x, 

I , A* ^ 1 + j^R^.a^a'' - ^VpR^^a'^a^a" , 

{IR -X)I- {^^d,R - \D'^F,^ - \D,X)a^I 

+ [—QD^DjyX + j^F^cyFj," — y^D ij.D" Fjjcr + 40 V^Vjy-R + ^20^ -^a''^ 

Tso^ ^^laup "T li^fjpT-Jrii, ZJrC^fjJrii, ) jcj (J 1 , 

at ^ iU-^R - Xf - \D^X + ^F,.F-^ 

+ -[^{Rfivcr pR^^'^ '' — RpuR^'^) + 30 V R)I , 

where a^{x^x') = g^^{x)d^a{x^x') and I{x,x') is the matrix giving parallel transport 
along the geodesic from x' to x defined by a^D^I = 0, I{x,x) — 1. F^^, is the usual 
field strength formed from A^. With the decomposition of Ga given in (1.5) extensive 
calculations at two and more loops have been undertaken using formulae for the singular 
parts of the products of G^ such as 

Similar relations involving derivatives may be found in ref. (15), in general at £ loops 
dimensional consistency requires a factor n~^^ . 






(1.7) 



In the neighbourhood of dAi the expansion (1.6) is no longer sufficient for reveahng 
all singular contributions in amplitudes at one or more loops. If a"(x, ic') is the geodetic 
interval on 9A^, as determined by the metric 7ij(x), then, for coordinates corresponding 
to the metric (1.1), assuming 

0(6^) , (1.9) 

and with the choice of gauge 
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Ga may now expanded as 



GA = i2 (^n + G^) + Gl , G^, G^ = 0(e-+2-'^) , (1.11) 

n=0 

with G^ regular both for x' ^> x and also on dM. for d = A. G^ correspond to the singular 
part of Ga 8iS x' ^ X while G^, which may be determined from the heat kernel expansion 
described in appendix A, are non singular when x' = x but are necessary in order to satisfy 
the boundary conditions and lead to additional divergences as the boundary is approached. 
If we define 

Gr{u)= ^^ ,, \ u^ + 2aY~^\ Gr,s{u)= dzz'Gr{z + u), ^^^^^ 

A-K-^"- Jo U-J-^j 

Gr,s+liu) +uGr,s{u) = ^sGr+l,s-l{u) + l6soGr+l{u) , 

then explicitly for n = 0, 1 we may obtain for general A and boundary conditions as in 

(1.2,3) 

G^ = Giiv) / , Gf = Go{v)uK,ja'a^ / , v = y-y' , 

G^ = Gi{u)V-I, P_ = 2P-1, u = y + y', 

Gf = Gi,o{u)2tlji + Go{u)uK,ja'a^V-i 

+ 2yy'{Go,o{u) K - 2G_i,o(tx) K,ja'a^)V-i 

+ (Gi,o(w) - Go,2{u))KV i + 2G-i,2{u)K,ja'a^ vi 

-Gi{u) 2Va'b,V I + Go,o{^) [Ay'VaW.V - Ay a'D.V V) I , 

where the coefficients, P, -i/;, Kij, are evaluated at x. /(x, x'), defined similarly to I{x,x') 
earlier, is the matrix corresponding to parallel transport along the geodesic in dM. from x' 
to X for the connection A = A|y=o and Di denotes the corresponding covariant derivative 
on the boundary, DiV = diV + [A^, V] . Kij is here the extrinsic curvature of the boundary, 
with coordinates given by (1.1) Kij = —^dy'jijly^o, K = ■j^^Kij. Clearly the presence 
of a curved boundary, so that Kij is non zero, is a significant complication. For a = 
the integrals in the definition of Gr,s in (1-12) are easily evaluated but in general the form 
given is sufficient for further calculations. To obtain (1.13) the terms Xln ^n in (1.11) 



(1.13) 



may be derived by expanding the singular contributions in (1.5) using (1.6) and (1.7) 
and the remaining terms Yin ^n '^^^ then be verified to be in accord with the essential 
Green function equation (1.4) and the boundary conditions by expanding about y = 0, 
A = Ao + Ai + . . . , An = 0(e''"2), regarding dy, Di as 0(e"^) in addition to y = 0(e) 
as in (1.9), and using results such as 

{u)u , 

(1.14) 
-d+1 , 

where Aq = —D^ — d^ and ViO"-^ = 5/ — \Rku^ a'^ a^ + . . ., for Riju the Riemann tensor 
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formed from the induced metric on dA4. The complete expression for the n = 2 terms in 
(1.11) may be found from the results obtained elsewhere for the heat kernel to this order, 
in appendix A, in (A. 10) we give the results for n = 2 involving X and t/j. These involve 
functions such as G2{u) where it is necessary to subtract a pole term at £ = as in (1.6). 

Instead of the result given by (1.13) an alternate equivalent form for Gf is given by 



(1.15) 



Gf = Go{u)u(l + ^)Kija'a^V-i 

+ (2yy'Go,o{u)V- + -^Gi,o{u)V - Go,2{u)v) [k - {d - 1)^^^ 

+ Gi,o{u){2iP +'^0^Kv)i 

-Gi{u) 2Va'D,V i + Go,o{u) {Ay'Va'D{P - Ay a'D^V V) I . 

This expression is in accord with the result that for a sphere of radius a when Kij = ^ij/a 
and also in the Neumann case for ijj — —^{d — 2)/aitis possible to find the Green function 
for —d^ in terms of the elementary flat space result for no boundary by the method of 
images. In this particular case (1.15) is independent of the special functions Gr,s- 

In any amplitude with superficial degree of divergence D then expanding Ga associ- 
ated with any internal line reduces D to D — 2n for the contribution of the Gn a^ term in 
(1.5) and to D — n for the terms G^, G^ in (1.11). For vacuum diagrams for renormalis- 
able theories in four dimensions D = A while for divergences associated with the boundary 
D-3. 

At one loop amplitudes are determined by the functional determinant of A and 
standard methods give 

-IndetA- / — Tr(e-"^), (1.16) 

Jo ^ 

and the singular part for d ^ A is then, using previous results [13,16] for the asymptotic 
form of Tr(e~'^'^) as r ^ 0, 

-IndetA--:^! / dvtY{\{\R- Xf + ^F^^F^^ + ^{R^^^pR^^'^'^P - R^^R^) 

e iDTT [ Jj^ 
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+ / dStvf^V-{ldr,R-dnX) + {2i;+lK)QR-X) 

JdM ^ 

+ J-KB +—K''^B- ■ -K'^^B- + -V D V F"^ 

XO ''Tin T^ 30 -'■''injn 90 *J 3 — -'-^*' -"■ n 

+ ^{{2V - f)K^ + (f P + 2)KK,,K^^ + (f - IV)K,,K^'^K,^) 
- ^D.VD'V K + ^D.VDjV K'^ - lD,VD'Vij^\ , (1.17) 

with dv = d'^x^ and dS = d'^~^5^^/^. 

Beyond one loop calculations become more involved. Using the representation out- 
lined above for the singular part of Ga in the neighbourhood of a boundary we have 
undertaken calculations at two loops for the simplest renormalisable field theory in four 
dimensions, i.e. purely scalar field theory with a multi-component scalar field (pi repre- 
sented by an action 



S{cP,g) 



[ dv i^d^cp.d^cp, + V{(P)) , V{(P) = i^9^Jke(p^<pJ(pk(pt + 0((/>3) , (1.18) 
Jm 



with g denoting generically the complete set of couplings, including mass terms, on which 
V depends, for vi{(/)) a basis of polynomials of degree 4 we may write V{(/)) = J2i 9i^i{4')- 
In order to take account the effects of the boundary dj^ we also add to 5" an extra surface 
term S{(f), g) which is a integral over the boundary of a local scalar formed from (j) and its 
normal derivative dnCJ) and which depends on additional couplings, or sources, g. We also 
follow the usual background field technique expanding 

(/> = V? + /, (1.19) 

for some fixed background field </? and assume boundary conditions on ip so that S+S has no 
terms linear in / when tp satisfies the classical equation of motion —'W'^ip-\-V'{p) = 0. The 
same boundary conditions are then imposed on (j), which therefore entails corresponding 
boundary conditions for the fluctuations /. The background fleld method allows an efficient 
procedure for calculating divergent parts of vacuum self energy diagrams as local dimension 
4 functions of </? and hence, since (/> and (p satisfy the same boundary conditions, obtaining 
the necessary counterterms Sc.t.{4>id)i Sc.tX4>^d^d) to ensure that the functional integral 

j rf[(/>] e-^"(<^'«'^) , ^0 = l^-' {S + ^c.t. +S + ^c.t.) , (1.20) 

defines a finite measure over fields (j) satisfying the required boundary conditions. The 
factor //~^ in Sq ensures that the couplings g, g, as well as the fields (p, retain their canonical 
dimensions although the quantum field theory is formally extended to d dimensions. 

In the next section we apply these methods to the situation where 5" is chosen so that 
't'ldM = 'P foi^ '^i(^) an arbitrary smooth boundary field. This corresponds to requiring that 



the quantum field / in (1.19) satisfies Dirichlet boundary conditions on dAi. In this case 
it is essential, besides the usual counterterms on M. necessary for finiteness, to introduce 
additional local counterterms on dAi linear in dn(p and also in K which are constrained by 
power counting to have overall dimension 3. These are calculated to two loops making use 
of the expansion (1.11) for the Green function Ga- The renormalisation group equations 
which now include terms involving the boundary operator dn(p are also discussed. Similar 
calculations are also undertaken in section 3 for the case when S is chosen such that dn4> 
is related to 4> on dAi and / satisfies generalised Neumann boundary conditions. In both 
cases the final results are also specialised to the 0{n) symmetric case, where there is a 
single coupling g and some results were also obtained by Diehl and Diettrich [5] . As usual 
with dimensional regularisation the renormalisation group equations determine the higher 
order poles in e of the counterterms in terms of the simple poles and these conditions are 
verified by our two loop results. The details of the calculation, for either the Dirichlet 
or Neumann case, are relegated to appendix B. They depend on the detailed form of the 
expansion given by (1.11), (1.12) and (1.13) or (1.15). A brief summary of the important 
results for the heat kernel for a general second order operator A of the form (1.2) with 
boundary conditions (1.3) in our treatment is given in appendix A. In section 4 consistency 
conditions are derived under local Weyl rescalings of the metric which in effect determine 
the K dependent counterterms in terms of those necessary, on fiat space and for a plane 
boundary, to define local composite operators on dAi. These results allow the derivation 
of a local renormalisation group equation which is discussed for the 0{n) symmetric field 
theory at the critical point where P{g) = 0. 

In section 5 we return to the situation for a general scalar field theory where the 
quantum field (p attains a fixed boundary value ip on dAt and the functional integral 
(1.20) defines a wave functional ^((^). Following earlier work by Symanzik [10], and in two 
dimensions by Liischer et al [11], this is shown, to two loops in a semiclassical expansion 
with the background field </? in (1.19) solving the classical equations of motion with the 
boundary condition tp = (p on dAt, to satisfy a local Schrodinger equation of the form 



^^^+^o(x)jvl>(^) = 0. (1.21) 

6t{x.) represents an arbitrary local deformation of the boundary dAt along the normal 
n^(x) and 7io(x) is a local Hamiltonian operator. To leading order in the semiclassical 
expansion TYq -^ ^ where the functional differential operator is 

'H = -\^ + V{^) + ir'^dpifid^^^ , (1.22) 

although Tio contains counterterms reflecting those occurring in 5*0 and also since the 
second functional derivative 5^/5(^(x)5(^(x') is ill defined in general as x' -^ x. Such 
divergences are demonstrated here to be fully resolved by use of dimensional regularisation. 
Furthermore the structure of TYq is shown to be determined in terms of Sq to all orders 
by use of the renormalisation group. Some mathematical details necessary to verify (1.21) 
are contained in appendix C. Finally in section 6 a few concluding remarks are made. 



2 Quantum Scalar Field Theory with Prescribed Boundary Values 

In this section we consider the quantum field theory defined by the action (1.18) with 
the additional surface contribution 



5((/>,(^)= / dS {<P - ip),dr,(i), , (2.1) 

JdM 

where (^(x) may be regarded as a source defining the surface operator dn4>{^) = dy(j){x)\y=Q. 

When (f) is expanded about a background </? as in (1.19) we may write S + S = 
Sc + Si + S2 + ■ ■ ■ with Sn = 0(/"') and Sc = S{lp, g) + S{(p, ip). It is then straightforward 
to see that 

Jm JdM 

SO that the appropriate boundary conditions to impose in a perturbative expansion are 

^\9M=^^ f\dM=^- (2-3) 

Also with this boundary condition 

S2 = \ I dv /,(A/), , A = -V^ + V"{^) , (2.4) 

JM 

SO that, acting on fields / satisfying (2.4), A is a symmetric operator of the necessary form 
(1.2), with X = V"{ifi), and the associated Green function Gij{x,x'), which corresponds 
to the propagator for internal lines in a perturbative expansion, is therefore defined by 
requiring Dirichlet boundary conditions. The higher order terms in the expansion in / 
as usual generate the vertices for Feynman diagrams in a perturbative expansion, in this 
theory they arise just from S'3, depending on V"'{(p), and 5*4 with no interactions restricted 
just to dM. 

For considering correlation functions of (f) in this quantum field theory with boundary 
conditions (2.3) the essential functional integral is 



The counterterms in 5*0 are chosen so that W{(p; J) is finite for arbitrary external sources 
Ji{x). For calculation it is convenient to define 

r(<^; ^) = W{^; J)- f dv Jm , (2.6) 

JM 

where ip is assumed to be an arbitrary smooth field on Ai, with boundary value ip, while 
J is constrained to be defined in terms of </?, with an invertible relationship J ^* </?. For a 
quantum field theory without any boundary it is conventional to choose J so that {(/)) = </?, 
to all orders in the loop expansion, and then r{(p) is the generating functional of connected 

9 



one particle irreducible amplitudes. This contains the only primitive divergences requiring 
renormalisation by appropriate counterterms. In a perturbative expansion F may be calcu- 
lated via the background field method, taking (j) = ip + f with (i[(^] = d[f], when to ensure 
(/) = it is necessary to require J = Jc + . . . for Jc(v') — A*~'^(~^^V' + ^'i.'fi))- However, 
in the presence of a boundary this condition on (/) necessitates that J{x) is singular as 
X approaches the boundary, assuming the background </? is smooth in the neighbourhood 
of d}A. This is a consequence of the fact that one particle reducible amplitudes also have 
primitive divergences, in the presence of a boundary, requiring corresponding contributions 
to Set.- For the purpose of the calculations of the counterterms required by the presence 
of a boundary undertaken here, it is simpler to choose J = Jc, without any higher order 
corrections, to define T{(p;ip) in (2.5) and (2.6). This choice, together with the boundary 
condition on (/?, cancels the linear term Si given by (2.2) in the expansion of the classical 
action although in this case it is necessary to include one particle irreducible graphs in the 
perturbative expansion beyond one loop. Since with dimensional regularisation there is a 
unique separation (</?) = {(fi)n.s. + (</')s.7 where (v')n.s.7 (</')s. are respectively non-singular, 
singular as the boundary is approached, an alternative possibility would be, to require 
(v)n.s. = but this is not pursued here. In this framework it straightforward to define a 
Scrodinger wave functional '^{(p) by ^((^) = expVF((^;0) = expr((^; (/?c) where (fc is the 
classical solution to — V^y^c + y'{.Vc) = prescribed by the boundary value (^. 

Following standard procedure the structure of possible counterterms is determined 
by power counting. As usual 

Sc.t.{(p,9)= I dv{\d^(t),A,^d^(t),+V,.iX(t))+C{X)) , (2.7) 

JM 

with Vc.i.{4>) a polynomial of degree 4 in (/> which is also taken to include terms linear in 
R which are quadratic in (f). C(A) contains all additional counterterms which involve no 
operators, or dependence on (/>, but are necessary for a curved space background, it may 
be written as 

C{\) = -{aF + bG + cR^) , A = (a, b, c) , 

r-n n^,p„ ^_^n «m- + (^_ 2)(rf- 1) ' 
G = R^^^^ R^j^pcr — AR^^ R^j^ -\- R , 

where A is determined perturbatively in terms of the dimensional coupling g. The addi- 
tional counterterms restricted to the boundary are of dimension 3 and have the general 
form 

S,.,X(t^,g,^)= [ dS{-ifc.t.^{0)^n<t>^ + p{0)K + C{X)) , (2.8) 

JdM 

using (pldM = 'f and also that with dimensional regularisation there are no linearly diver- 
gent terms which are 0((^^). Assuming minimal subtraction Set. and Set. contain only 
poles in e. C, with A denoting the coefficients in some basis, includes all terms which are 
proportional to any dimension 3 scalar, independent of ip, formed solely from the extrinsic 

10 



and intrinsic curvature on dAi, such as appear in the surface terms in (1.17) not involving 
t/j and X. 

Due to the complexity of the calculations necessary to determine such contributions 
we do not consider them here but focus primarily on ipc.t., which has in effect been calcu- 
lated previously to two loop order by Diehl and Dietrich, and also on p. Of course Set. in 
the form (2.7) has been calculated extensively before at two and higher loops although it is 
necessary to redetermine all contributions to A since any integration by parts of derivatives 
may lead to surface terms. It is also necessary to determine any divergences proportional 
to total derivatives which are usually discarded. Using the background field method 5'c.t. 
may be calculated with dn4> -^ dnV- It is easy to see that ipcxXv) is linear in (p and 
depends in general on V"'{ip) while p{ip) is quadratic in (p and is proportional to V"{ip) or 
V"'{(pY ^'^d no other ip dependent terms are required in (2.8). 

In general for a perturbative expansion of the functional integral (2.5) with a back- 
ground field (fi 

r(^;v^) = -5o(¥^) + $^rW(^;v'), (2.9) 

e.=i 

where V'^^^ip) is the connected amplitude at £ loops and 5'o((/?) contains the necessary 
counterterms to ensure that r((^; </?) is finite. At one loop 

r(i) =-ilndetA , (2.10) 

and using (1.17) it is easy to see that to ensure cancellation of the poles in e it is sufficient 
to take in (2.8) 

^'clM) = "T^ ^y^j^i^ ' P'''(^) = -T^ i^(^) • (2-11) 

Beyond one loop the calculational details are non trivial and are mostly relegated to 
appendix B. The essential two loop vacuum diagrams are shown in fig.l and accordingly we 
decompose F*^^) = F^ + F^ + . . .. Of course at two loops there are one loop sub divergences 
which are removed by one loop counterterms. For the scalar field theory defined by (1.18) 
and (2.1) these are given solely by 

V}1] - c^'^R' = ^ i V,,V,, , y,,(<^) = V:;{cl>) - \R5,, . (2.12) 

The amplitude for fig. la, with appropriate counterterms, is then 

Fa = - — / dvgijkeGij\Gke\ - ^ - / dv gijuVij{ip)Gu\ 

W / 4^- ~ " _ , _ , (2.13) 

= g / dvQijU L^g 2 )2 ^iM^k(.{^) - p''Gij\Gke] 

where in general G\{x) denotes the the coincident limit, x' — > x, of the Green function 
G{x,x'). In (2.13) we have used 
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where, for manifolds without boundary, Gij\ is finite for d = 4 from (1.5), (1.6) and (1.7). 
Similarly for fig. lb, with its counterterms, 

(2.15) 
Using (1.5) and (1.8) it is straightforward to see that 

r^ - 0^2-, IJdv y- (v.) (v^ + ii?) v:;^ (v.) 

^ -e \ . (2-16) 

-J^ 1(1-1, 



+ 7T^Z27V2 o (1 - ¥) I d^ V:jk(y^)V:;',{y.)V,,{^) 



and together with l/e^ term in (2.13) this determines the required two loop counterterms 
for 5'c.t. as usual. 

In addition it is essential to consider also the one particle reducible diagram fig.lc 
for which 

Tc = ^ j j dvdv'G,u\Vi;'^{^)Gu'Vir^,^,{^')G,,,,\ . (2.17) 

For no boundary T^ is finite while here it is necessary to allow for the extra one loop 
boundary counterterms given by (2.8) and (2.11) which give the one particle reducible 
diagrams in fig. Id and 



^'^ " 4 16^ j j dSdv'Vl!j'^{^)d^Gu'\y^^Vir^>A^')Gj'k'\ 



(2.18) 



When sub divergences are subtracted in one particle reducible amplitudes then according 
to the usual lore of renormalisation theory there is no remaining overall divergence. If 
Gii'{x,x') in (2.17) and (2.18) is separated into a singular and a regular part then by 
integrating by parts from (2.17) and (2.18) we find 

(Fc + r,)-s = ^ / / dvdv' D,,V{;',{^) G^fV^/^,,, {^')D,,,, , 

^ ^ ^ _^ (2.19) 

Djk = G,,\ + Y^^^jk{-S'{y) + \K5{y)) , 

where Djk{x), as in (C.20), is the dimensionally regularised form of Gjk{x^x) with poles 
in e representing both the usual local and also boundary divergences subtracted (the 5{y) 
term is irrelevant in this case assuming G^^F also obeys Dirichlet boundary conditions and 
vanishes at y = or y' = 0) and hence, so long as G^^^f^x.x') is a sufficiently smooth 
function on A^ x M., (2.19) is finite as £ — * 0. However, as pointed out by Symanzik 
[10], for the full Green function Gii>{x,x') there are additional divergences in (2.17) not 
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removed by the counterterms (2.18) and in general the manipulations leading to (2.19) 
are not justified. An additional complication in this case is that in (2.18) dnGu'dn has a 
non integrable singularity as x' ^ x which is not regularised by analytic continuation in d 
as usual with dimensional regularisation. Nevertheless elsewhere [17] we have shown that 
there is a well defined prescription for handling this which we follow in detailed calculations 
in appendix B. 

The results of appendix B allow the divergent terms as £ —> on the boundary to be 
readily calculated for each contribution. For Fa from (B.4,7), and also using (B.8,10), we 
get, keeping only terms relevant for determining ipc.t. and p in (2.8), 

Ta - (y£f i jdSg.,,,{{l - \e)d^Vl'^{^)l^^ - |(1 - le)Vl^{^ k] . (2.20) 

The amplitude F^ also has, apart from (2.16), a contribution to the divergence on the 
boundary given by, from (B.19) and (B.20a,..f), 

^ ' ^ ^ dvV.,V''iV^V''A 



(167r2)2£ 24 7 ^^ '^^'^Jk^ Ujk 

+ j£^ IJ^H^^- -e'^^^^(^^9r.v{;',{^)\^^^ - i(i - le)v:;^{^v!;',{^ k} . 

(2.21) 
For the remaining one particle reducible amplitudes we may use (B.33) to give 

(2.22) 
From these results we may then read off 

'^^''■' = ~ (167r%)4 i^^ " '^^^^^^Sjkli + 1(1 - \e)9^uyM " ^(1 - ^)9^nkyu^ , 

(2.23) 
With 5*0 defined as in (1.20) we may write 

Soic/)) = 5((/>o, go) + S{c/>o, ^o) + C{p) + C(A, A) , 

r ^ . r r . . (2.24) 

C{p) = ij-^ dSp{ip)K, C{X,X)=p-' dvC{X) + p-^ dSC{X), 

JdM Jm JdM 

in terms of the expressions for 5", 5" given by (1.18), (2.1). C includes all necessary (/>, (^ 
independent counterterms not calculated here. To obtain (2.24) it is necessary to require* 



p—^'Z,j(Pj, Z={l + A)-^ 



0Oi = P ^'^{Zij0j+Z ij0c.t.j) ■ 



* Our treatment at this point differs from Symanzik [10] and also Diefil and Dietrich [5], although 
results to two loops are equivalent with the latter. 
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The renormalisation group equations are derived as usual by requiring invariance of 
physical amplitudes under rescalings of the arbitrary mass /U. To obtain these it is sufficient 
to require 

^i^So = h- f dv{^(P)~)So = -C(/30 -C(/3\/3^) , 
da \ L. 6<hJ ^2.26) 



v = ^^ + ^y.^+ [ dS^^— 



dn dV Jqm S0i 

where /3^ ■ d/dV = P^jd/dgi and pp,P^,(3^^ are finite. By using (2.26) with (2.5), and 
discarding a total functional derivative in the functional integral, it is easy to find the 
renormalisation group equation, 

(V + J dv{^J)~)w{if;J)=C{n+CW\p'^) . (2.27) 

As is conventional for field theories without boundary in dimensional regularisation on the 
basis of minimal subtraction 

/3^((/>) = e{V{cP) - lv',{<p)<p,)+(3''i<p) , % = -i£5,, +7., , /3^- ^Z = Z^ . (2.28) 

The P functions associated with the boundary that may be obtained by the calculations 
here are (3'^ and pp which from (2.26) are determined by 

/3^- gy+^'^ji^)-Q^. - i^)^Orm = , /3^(^) = i£^ + /3n^) , (2.29a) 
/3^-^ + /3^,(^)^-£)p(^)= -P^m. (2.296) 

From (2.29a) we therefore find to two loops, since /9^^^) = ^V-jV-j/{16n'^), 

L 1 

167r^ 

+ riQ^2)2 Y^{^^ijki'f)9jkU + '^dijklVjkii^) - ^9ijjkVka{ip) - 9ik£m9jk£m'fjJ ' 

(2.30) 

using (2.11) and (2.23) with the two loop result for Z^- = —^gik£m9jk£m/{^Q'^'^)'^£- As 

usual with dimensional regularisation the double poles in (^^ ( are determined by (2.29a) 
which provides a significant consistency check on our calculations. In a similar fashion 



(^^^i^) = -T7r^lKji^) 



+ 7T7rz2V2 i^\^9ijkkyij{^) + ^yijk{^)yijk{^) - yijM)yikk{^)\ ■ 
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For specific application of tfiese results we focus on the 0{n) symmetric potential 
V{(j)) = ■^g{4>^Y + \m'^4)'^ : for (pi an n component field. In this case we may write 

^%ip) =fi0,, I3P{0) = Ip0^ + Pmm' , (2.32) 

so that t) is the anomalous dimension of the boundary operator dnCJ). From (2.30) and 
(2.31) to two loop order 

f}= -\{n + 2)u + ^(n + 2)^2 , 

1 n (2.33) 

p= -i(?i + 2)w + ^(?i + 2)(n + 5)w2, p^ = ^^-(-l + |(?i + 2)w) , 

for u = (//(IGtt^). For comparison with Diehl and Dietrich [5] f) = ^rji + 7 where for this 

model 7ij = 'y6ij and 7*-^-' = -^{n + 2)vF'. 

As an illustration of the renormalisation group equation we apply (2.27) in this 0{n) 
model to the expansion of the field operator (t){x) in the vicinity of the boundary. Defining 

then we may write, for y ^ and neglecting to, 

((/.(x, y))j^ Coipy) <^(x) + Ciipy) y{dn(t>{^))j + C2{py) 2/i^(x)<^(x) . (2.35) 

From (2.27) the leading singular parts of the coefficient functions satisfy 

[p— + i3{g) — + -f{g) + ?)((7)) Coipy) = , 

(^^ + /3(^) ^ + 7(^) - vig)) Ci (py) = , (2.36) 

[pg- + kg)g- + ^(9) + v{g))c2{py)= - p{9)Ciipy) , 

where manifestly for g = Cq = Ci = 1, C2 = 0. This result defines an expansion for the 
operator (j){x), analogous to the usual operator product expansion, in the neighbourhood 
of the boundary. At the infra-red fixed point corresponding to (3{g^) =0 

Coipy) oc (^y)-^*-'^* , Ciipy) oc (^2/)"^*+'^* , C2{py) ~ -^Ci{py)+0{{py)-'^*-^*) . 

(2.37) 
From (2.33), and using that with minimal subtraction the critical coupling is given by 

3 9 

«* = 7 ^ e + ^T:T(3n + 14) e^ + 0(£^) , 

(?i + 8) (n + 8)3^ ' ^ ' ' 

then the critical indices are as usual expressed as an expansion in e 



1 (n + 2) ^ ^ 1 n + 2 1 (n + 2) , „ „ , . 

7^ = - ; ,; e^ , TO = ^ £ - - ; ,; 17n + 76 e' 

' 4 (n + 8)2 ' '* 2 n + 8 4 (n + 8)3^ ' 

p* _ 2 _ 1 n + 4 

17* 3 3 n + 8 
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(2.38) 



3 Quantum Scalar Field Theory with Free Boundary Values 

In this section the essential quantum field theory given by the action (1.18) is ex- 
tended by the surface term 



S{(t^r9)= / dSQ{(t^), (3.1) 

JdM 

where Q{(p), depending on couplings g, is a scalar and there is no dependence on the normal 
derivative dn4' as in (2.1). For renormalisability Q{(p) is at most cubic in (p although if 
(J) — * —(/) symmetry is imposed at leading order Q{4') — 0((/)^) and this is stable under 
renormalisation. 

When (j) is expanded about a background (p then to first order 

^1 = / dvf,{-V^^, + V',{y^)) + [ dSf,{-dny^ + Q'{^)), , (3.2) 

Jm JdM 

SO that the appropriate boundary conditions, assuming that they apply to (f) as required 
for invariance under the shift symmetry generated hy 6ip = —df, which should be imposed 
are now 

(a^<^-Q'(</>))4_^=0. (3.3) 

To second order in / the action is again of the form (2.4) where A is now a symmetric 
operator acting on fields obeying the generalised Neumann boundary condition 

{^nf^-Q'lJ{^)fJ)\QJ^=0. (3.4) 

If Q((/>) = \gijk4>i4>j4>k + . . . there is an additional surface interaction in this case 

"^3 = -777 / dSgijkfifjfk , (3.5) 

-■-^ JdM 

arising from (3.1) and also 0(/^) terms in the expansion of (3.3) but we do not undertake 
any calculations involving this term here. 

The basic functional integral remains as in (2.5) but now integrating over fields 
satisfying the boundary condition (3.3) and defining now W{J). The necessary local 
counterterms on TW for finiteness of W are just as usual for the scalar field theory without 
boundary, as in (2.7), but on the boundary dA4 they may now be written as 

^c.t. (</>, 9.9)= [ dS (Qct. {<P) + C'( A)) , (3.6) 

JdM 

where C{X) represents contributions independent of (p which include all terms of the same 
form as those present in C for the Dirichlet case in (2.8) but also, for gijk = 0, any terms 
involving dimension two scalars formed from the metric, such as Rnn, KijK"^^ or K^, and 
proportional to Q" . Qc.t.{4') then contains the remaining counterterms independent of the 
metric and also those linear in K. In the background field method Qc.tXv^) ^^^J be calcu- 
lated directly and involves, for c/ijk = 0, just terms with the general structure V"'{ip)Q'{Lp), 
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V"{f)Q", V"'{^pfQ", Q"^ while terms proportional to K may contain V"{f), V"'{f)'^ or 
Q"^. Clearly Qc.tX'^) remains 0((/7^) in this case. 

The one loop amplitude is still determined by the functional determinant of A as in 
(2.9), although now with the alternative boundary conditions (3.4), and (1.17) gives, after 
using (3.3) to eliminate c^nV' terms, 

Qi'l = Y^ {-kv^ljQ'j + nv"Q"h - liQ"'h + liQ"'h K - \vll K) . (3.7) 

The corresponding result for C^^' is easy to read off from (1.17). 

At two loops the essential amplitudes are still given by (2.12), (2.14) and (2.16), 
if dijk = 0, although Gij{x,x') is now the Green function for A subject to boundary 
conditions given by (3.4). The one particle reducible diagrams in fig. Id, with the one loop 
boundary counterterms determined by (3.7), now give 

r. = ^ 1 1 dSdv'Q^^lMGu'\y=o^l''M^')Grk'\ 

-err (3-^) 

In a similar fashion to (2.19), if Gii'{x,x') — > G^^F(x, x') which is assumed to satisfy the 
same boundary conditions, then combining (2.17) and (3.8) gives the finite expression 



(rc + r,)-s = ^ j j dvdv'D,uVi;',{^)G'^fVir^,,,{^')D,,k, 



n 



-e 



(3.9) 



^^■^ = ^^-^1 + ikr. i^o^ ^'^y^ + (4g;', - 5,4^) 5{y)) 



IQ-k'^s 



since Djk is in accord with (C.20) for this case as well. Since Qc.lX4') is quadratic in (f) 
there is an additional diagram fig.le which corresponds to 

Te = — j dSQ^^ -J Gij\ 

^ ~ 4 T67r% / '^^^^^^^ i^kiS'iy) + {AQ'^e - 6k£^K) 5{y))Gij\ . 

The results of appendix B again allow the divergent contributions on the boundary 
to be calculated for terms independent of and linear in K. From (B.4,6) we get 



Ta + Te ~ ^yI^ JdS {g.jkk{2{Q'%j - (1 - e){Q"V") 



ij 



- i(2 + e)Vl^^[ - 1(1 - e){Q'%, K-\{1 + le)Vl', K^ 

+ 9^Ju{{^ - e)Q'I^Vl', - AQ'l^{Q"^)u + |(2 - \e)Q':^Q'U k) } . 

(3.11) 
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Also from (B.19) and (B.20a,..f) 



Th ~ - ...^... ^ dvV„V''iV^V'" 



(3.13) 



+ (y£)2 jds{-\{l- eWl^, g,^,, Q[ + {l-\e- \'k\)VZ V^'^, Q'l, (3-12) 

For the remaining one particle reducible amplitudes we may use (B. 27, 28, 31) to give 
Tc + Td ~ ^g^2g-)2 / "^-^ 1 1(1 + ^y^iik 9jjke Q'i + k^iik ^j'j'e Q'U 

~ 2^iU ^jki Qjk + 12 (-*■ ^" 6^)^iik ^jjk ^ 

Combining these we therefore obtain the somewhat lengthy result 

QS. = (Y^^ {9^M'^{Q"%J - (1 - e){Q"V'% 

- 1(2 + e)Vl-^', - 1(1 - e){Q"\, K-\{1 + leWi; K) 
+ g^Mi^^ - e)Q':^Vl', - 4Q'r{Q"')ki + |(2 - ^e)Q'^^Q'^, K) 
+ g (1 + £)V^ik 9jjke Qe + s^iik ^jje Qki ~ 'qYui ^jke Qjk 

+ T2 (-*- + e^Wiik Vjjk -^ ~ 2 (-*- ~ ^Wijk 9ijke Qi 

+ (1 - i^ - h"^)v:;.'i ^^^ Q'!j - 1(1 - 1^^ - h"^)v:;^ v^jk k} . 

(3.14) 
To obtain the renormalisation group equations we define, as in (2.24), 

5o(</>) = S{(po, go) + ^(</>o, ^o) + C(A, A) , 
Qo(</>o)=^-^(Q(</>)+gc.t.(</>)) , 

for C(A, A) formed from the counterterms C(A) and C(\) in this case. Since Qc.t. contains 
terms involving K, which is in general x dependent, it is necessary to regard the couplings 
^ in Q as also x dependent to ensure multiplicative renormalisability or that Qq is obtained 
from <5 by ^ -^ ^o- In this case instead of (2.26) we may then define the (3 functions by 



7 /* X 



(3.16) 



In particular f]'^ is determined by 

(/3^- ^ +/3«- ^ - {7<Ph^^-e)Qo = , (3Q{<P)= e{Q{cP) - iQ ',(</>)</>,) +/3^(</>) .(3.17) 
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Although Qo? 13^ are both given as hnear in K it is possible to generate K^ terms in (3.17), 
with the assumed structure of counterterms even if Q'" = 0, but these may be neglected 
since they are not calculated here. To two loops the P function becomes 

/3« - q:(7</>). = ^ {-^v:i;q^ + 2{v"Q'% - i{Q'%, + mQ"\^ - v:!)k} 

"*" 7Tr~~2Y2 \9ijkk['^\Q ^ )ij ~ j^ijeQe) ~ '^dijkiQij^ke 
+ 4^iik 9jjke Qe + ^ijk 9ijk£ Q^ ^ (1 + 3"^ J^iki ^jki Qij 

+ gdijkk (6(<5 )ij — ^ij)K ~ gdijkeQijQke K 

+ iei^lLV^jk + (^+^^')V:;'kV^;^K] . (3.18) 

Instead of (2.27) the renormalisation group equation is now 

{V + j dv{^J),^)w{J)=CiP\p'^) , (3.19) 

with V as in (3.16). 

For the 0{n) symmetric case considered at the end of the last section, where there 
is just a single coupling g, we may write 

Q{cP) = ci(/>2 - h,<p, , (3Q{<p) = (3,\cj? -f,k<p, + X , 20) 

(3c = r]cC + pK , x = ^mcm^ + acC^ + {pmrri^ + PcC^)K + u dnrri^ , 

where the term proportional to dnW? is present if w? is dependent on x and is discussed 
further in the next section. The renormalisation group equation (3.19), setting the curva- 
ture on M. to zero, becomes 

= p-^ f dv^pm^+fj.-' f dS{x + C{(3^^)) , (3.21) 

Jm JdM 

where P'^ {(p) = ^P(p'^ + \^mm'^4>^ + ^pm'^. 

The results of our calculations then imply to two loop order 

rjc = lin + 2)u + ^{1 - 47r2)(n + 2)u^ , ?) = -i(n + 2)u + i|(n + 2)u^ , 

In/ 2/ \ \ 

P"'^~m^3V^^^''^ ^V ' (3-22) 

^ \ . 1 . 1 4n 

where r]c and f) are in accord with the results of Diehl and Dietrich [5] with tjc ^ Vc, V = 
^??i + 7- In a general regularisation scheme c is linearly divergent but as usual such 
contributions are absent with dimensional regularisation. 
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4 Conformal Invariance and Consistency Relations 

The essential new results of the previous two sections has been to calculate, to two 
loops, those counterterms necessary for a renormalisable quantum scalar field theory which 
depend on the extrinsic curvature K of the boundary. In this section we show how consis- 
tency conditions resulting from considering arbitrary local Weyl rescalings of the metric on 
M. and hence also of the induced metric on dM. may determine the K dependent countert- 
erms without undertaking any calculations beyond those for a fiat spatial background and 
with a non curved boundary. The method is an extension of that used to derive consistency 
relations which determine the dependence of counterterms on the Riemann curvature for 
general renormalisable field theories on manifolds without a boundary [18]. 

Under an infinitesimal Weyl rescaling the essential results here are 

d^g^"" = 2ag^''', S^R = 2aR+2{d-l)V^a , S^n^ = an^ , 6^K = aK+{d-l)d^a . (4.1) 

The consistency conditions are derived by considering relations involving the finite local 
scalar composite operators constructed from the basic fields. To define these the essential 
couplings (77, gj parameterising the theory are extended to arbitrary local functions gi{x), 
as well as gi{'x.), or effectively V -^ V{(/}, x), Q ^> Q{.4>^ ^), so that local operators may be 
obtained by functional differentiation. The quantum field theory remains finite, in renor- 
malised perturbation theory, so long as appropriate additional counterterms depending 
on derivatives of the couplings are introduced. In particular on dM. there are possible 
counterterms depending on d^gj. 

Here we consider initially the basic scalar quantum field theory with classical action 
(1.18) with Neumann boundary conditions so that there is also a surface term as in (3.1). 
If V{(f)) = J2i 9i'^i{4') and (5((/>) = X]/ QiQii'P)^ fo^^ some basis of quartic polynomials vi{(j)) 
and quadratic polynomials qi{(j)), then we define Agj in terms of AV{(j)) = 4V {(j)) — V'i{(j))(l)i 
and Agi by AQ{(j)) = 3Q{(j)) — Q'i{(j))(j)i so that if Agj = Ajgj and Agj = Ajgj then A/, 
Aj are the dimensions of the associated couplings. Defining 



A^= [ dva(2g^''-^+Agi^) , A^ = / dSaAgi^, Af = [ dva<P~, 

(4.2) 
then for d = 4hy simple scale invariance it is easy to see from (1.18) and (3.1) that 

{A^ + A^ + Ai){S + S) = -^ f dv(j)^V^a-l f dS cj)'' d^a . (4.3) 

Jm JdM 

The V^o", dn(T terms appearing on the r.h.s. of (4.3) may be cancelled by introducing 
appropriate contributions involving R(p'^ in S and Kcf)^ in S which ensures a classically Weyl 
invariant theory under local rescalings of the metric as in (4.1) if there are no couplings 
with dimension, Agj = Agj = 0. Later we show how this may be extended to the quantum 
theory at a critical point (3{g) = 0. 

In the full quantum field theory, for general d, where S + S ^' Sq the corresponding 
equations representing local changes in scale involve the f3 functions of the various couplings 

20 



and the anomalous dimension of (j). Defining in addition to (4.2) 

Af = / dva(3f^, Af = / rf5a/3f-^, A^^ = / dva{^<p)~. (4.4) 

Jm (^91 JdM ^91 Jm ^(Pi 

then the local renormalisation equation may be expanded in the form 

(A^ + A^ + A^ + A2^-Af-Af)5o = ^-^/ dvaC{(3^)+^-' [ dSaC{(3^) 

Jm JdM 

~fi-^ f dv{d^aZ^' + \/^aT)-n-' f dSdnay, (4.5) 
Jm JdM 

where Z^, T and y may depend locally on (j) but not its derivatives. Since gi is required 
to be X dependent for a consistent treatment it is necessary to include additional terms in 
(3^ {(/)) = X]/ (3i''^i{(p)i ciiid also in C{(3^), depending on d^gj which have been discussed in 
detail in ref. [18]. To absorb all counterterms by a redefinition of couplings in this case 
it is also necessary to introduce an external gauge field coupled to (j) but this is omitted 
here since it is unnecessary for obtaining the essential results of our discussion. When a 
is constant then (4.3) is equivalent to combining the usual renormalisation group equation 
(3.16) describing variations in ^ together with a simple scaling relation. Since Sq defines 
a finite theory for arbitrary g^u{x) and gi{x), 5'/(x) the l.h.s. of (4.3) represents a linear 
combination of finite local operators and hence the new quantities Z^, T and on the 
boundary y appearing on the r.h.s. of (4.3) must therefore also be expressible in terms of 
finite local operators. In the first paper of [18] Z^ and T have been calculated to two loops 
and consistency relations following from the requirement of finiteness have been analysed. 

Here we consider just the surface contributions depending on dna. To calculate y in 
(4.5) it is convenient to decompose the counterterms for Q in the form 

Q,,,,{4>) = q{4>) + p{4>)K + wM)d^gi , (4.6) 

where the last term is new in this case when the couplings are allowed to be x dependent. 
It is then easy to see from the defining equations (4.3) and (4.1) that in general 

y{ci>) = u{<p) - (3 - e)p{cp) - wM){/\gi - /3f ) , (4.7) 

if u{(f)) is the surface contribution resulting when the volume integrals are integrated by 
parts to achieve the desired form (4.5). 

The results for p to two loops are contained in (3.7) and (3.14). For uJi{(j)) then a 
straightforward extension of the calculations in section 3 gives 



wY^d^gi = - ^^-^^ \d'ryll , (4.8a) 



^f^^r.91 = (16^2^)2 1-1(1 + ^^)9^Jkk^'nV,'; + §(1 + sWZO'ry^^k 

- i(l - e)V^;',d'ryi;^ - Id^g^.^k V^/; + (l + e)d^g,,kk{Q"\j] . (4.86) 
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where d'n = n^d'^ with d'^ denoting the derivative at constant (j). To obtain w((^) then 
taking into account (4.3) and the derivative terms which appear at two loops in the volume 
integral as in (3.12) to this order gives 

«(</>) = i(l - \e)ct? - ^Y^^ ^(1 + \e) Vl'^MW^^^ + • • ■ ■ (4.9) 

As already remarked the consistency conditions stem from the fact that each term 
in (4.5) must separately be finite as operators, or on insertion in the regularised functional 
integral defined by 5*0. Since arbitrary finite local composite operators corresponding to 
/i((/)) = hjqj{(/)) on dA4, for h{(/)) at most quadratic in cj), may be defined by 

hijJpr^So, (4.10) 

then we may write 

/i-^W)=/i7T^5o, (4.11) 

ogi 

for suitable h{(p), non singular as e ^ 0, in this case. As a consequence of minimal 
subtraction (4.11) all terms involving poles in e in y{(p) are determined in terms of h{(p) 
and the Q dependent counterterms in Sq. h{(/)) also is determined by the 0(1) and 0(e) 
terms in (4.11). It is easy to read off from (4.7) and (4.9) the lowest order contributions 

h{<p) = i(i - ie)<p' + ^ I iv::{<f>) - Q'l.imiM)) + • • ■ ' (4-12) 

where the one loop contribution arises from ep^^\(j)) in (4.7) and also from the —^ecf)^ term 
in h{(t)) in conjunction with q^^\ (4.11) then provides the crucial consistency relation which 
in effect determines p in (4.6) in terms of lower order results in the loop expansion. We 
have verified that this confirms the one and two loop results of section 3. As an example 
from (3.7) 



^p'"' = -^MQ"")- - ^n ^ ^''' = T7^A^iy"Q"h - \v::^Q', - W)^) 



and from (4.8a) 



^f'^9i = -^iyii-\vi:'^c^^ 



so that it is easy to verify (4.11) at the lowest non trivial order since 

?,p^'^ + w^P^g: = -h',^q^'\ (4.13) 

ogi 

for h^[(j)) = i(/)^. The same procedure has also been used to check p^"^' using h in (4.12) 
although many more terms contribute. 
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From (4.5) we may derive a local renormalisation group equation for the functional 
W. Restricting to the 0{n) symmetric model considered at the end of sections 2 and 3 and 
also assuming the coupling g = g^ is at the critical point where /3((7^) = then defining 

A„ = j^i.a (2,,-^ + (2 - -,^.W^ + (3 - . - i.) J. A) 

/* c r 

+ / dSa({l-7]c,)c— + {2-^e-fi,)h^^, (4.14) 

JdM ^ "C Shi^ 

A;= / dv{V^a+^aR'-f^,)T—^+ dS {d^^a + a K'r]c.)9— , 
Jm ^^ JdM ^c 



where 

the local renormalisation group equation has the general form 

(A^ + A'^)W{m^, c, J, h) = X{m^ - \tR\ c - OK') , (4.15) 

with X a volume or surface integral of local scalar functions of m^ and c proportional to 
a or its derivatives of dimension 4 or 3 respectively. The arguments of X are written in 
the form shown in (4.15), without any loss of generality, for later convenience. From (4.5) 
and (4.10) the coefficients 9, ic^ ^m are determined by writing h{(/)) as 

h{(p) = ^e (p"^ + LcC'^ + Lmm^ , (4.16) 

with the coupling g evaluated at the critical point. Using (3.20,21) with g = g^ and the 
consistency relations (4.7), together with previous results in the case of theories without 
boundary [18], gives 

X(m2, c) = - ^i-' f dv {ai^pm"^ - (3aF - (3bG) + (VV - ^aR'{-f^, - e))km^] 
Jm *- -" 

- H~^ dS la{acC^ + a^cm'^ + PmK'm'^ +udnm'^) (4-17) 

JdM *- 

- dn(J imrri^ - {dn'y - crK'{2r]c^ - e))LcC^ + . . . j , 
with 

P'm = (d- l)pm + Oam = {im* " £)''m " {2 - 'Jm*)'^ ■ (4.18) 

The neglected terms in the surface integral in (4.17) are linear in c and involve invari- 
ants constructed from Kij and the Riemann tensor on Ai for which the corresponding 
counterterms are not calculated here. 

The restricted form of the functional operator A'^r in (4.14) and also X in (4.17) 
is a consequence of the consistency relations such as given by (4.7) and (4.11) at the 
critical point. Alternatively the expressions given by (4.14) and (4.17) may be regarded 
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as determined by the integrability condition [A^^ + A'^r, A^r' + A'cr/]VF = 0, without any 
commitment to a specific regularisation scheme [18]. Such conditions also require the 
absence of a R^ contribution in addition to the F,G curvature terms in (4.14). From 
(4.12) and (4.16) we may find to one loop 

= l-'^e + l{n + 2)u., ^c = -^ j^ , ^- = ^ J^ ' (4-19) 

while from (4.8a,b) to two loops 

1 TL 



UJ 



In terms of the definition of /3c in (3.20) then (4.14) requires that at the critical point 
{d — l)p — —TjcO which along with (4.18) may be verified with the results (3.22) and 
(4.19,20) using 7^,* = \{n + 2)^* + 0{u1). Even at the critical point W is arbitrary up to 
the addition of local functionals of m? and c. For 



^i-' f dS {rK'm^ + sdr,m^) - p-^ [ dvtR'r 

JdM Jm 



then 

dLrn=r+ {2-^rn*)s , SuJ = {^rn* - £)s , Sk = 2t , 

SO that in the e expansion it is possible to set irn and k to zero but not cD. 

Given the form of the renormalisation group equation (4.14,15) at the critical point 
it is straightforward to see that it may be rewritten as 

A^W{m'^ + \tR', c + eK\ J, h) = X{m^, c) . (4.21) 

Manifestly W{^tR' ^ 9K' ^ 0, 0) is invariant under local Weyl transformations, apart from the 
contributions involving X on the r.h.s. of (4.21). This is a consequence of the result that, 
at least for simple field theories, invariance with respect to constant scale transformations 
leads to symmetry under the full conformal group [19]. Using previous results [18] in the e 
expansion r = 1 — \e+^{n + 2)u\ + 0{e^) so there is a modification of the usual coefficient 
of the 4>^R term in a free scalar theory for conformal invariance in the interacting case at 
three loops. 

From (4.21) it is straightforward to determine the behaviour of the composite oper- 
ator (j)^ in the neighbourhood of the boundary at the critical point. Defining 

^^-^{ct^\x)) = -J^^W{m' + \rR'M',Q,Q)\^^^^ , (4.22) 

then from (4.21) we may obtain, with the standard coordinates x = (x, y), 

(I dva2g^''^-Xma{^,y)){<P'{^,y)) 

^ -co a°(x) {d'{y) - {xm - l)i^'(x) d{y)) ^^'^S) 

- im{dncr{^) - {xm " 2) a^ {^) K ' {^)) 5 (y) , 

24 



where Xm = d — 2 + '^^^^ and (J^i^x) = (t(x, 0), 9^(t(x) = dya{'s.^ y)\y=o- Using that under 
a Weyl rescahng as in (4.1) 

d^y = -aO(x) y - \dM^) v" " \^l<y{^) ?/' + ••• (4-24) 

then we may solve (4.23) for the leading singular behaviour as y ^ 0, 

((/>2(x,y)) ~ [y— + \x^K\^)y—^^)C 






V Z Xrn. ^ 



For y > this is in agreement with a general formula of Cardy [9], apart from a sign. 
The S'{y), d{y) terms reflect the regularisation of the singularity for Xjn -^ 2, when it is 
necessary that C -^ —u. 

A similar discussion to the above case can also be essentially carried out in the 
Dirichlet case, corresponding to the calculations of section 2, with similar consequences 
but it is not presented because of some technical problems. 
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5 Schrodinger Equation 

The results of section 2 verify the feasibihty of defining in perturbation theory a finite 
wave functional "^{(p) for renormalisable scalar quantum field theories on a manifold Ai 
where the quantum fields (p are specified by an arbitrary smooth field ip on the boundary 
dAi. For such wave functionals there is a corresponding Schrodinger equation which 
describes the variation of \1/ under smooth local deformations of dAi. A discussion related 
to that presented here was also given by Symanzik [10] although only constant global shifts 
of the boundary were considered. The relevant Hamiltonian is a functional differential 
operator requiring careful regularisation as is discussed subsequently. 

The wave functional ^ is defined formally by 

ii(0) = I d[(t)] e-^^o^-^) , (5.1) 

where we have introduced Planck's constant % in order to set up a semiclassical expansion 
in powers of %. In general the boundary of A1 is supposed to be determined by a;^(x) for 
X coordinates on 9A^ and we then consider variations along the normal to 9A1 by taking 

5x^(x) = -?i^(x)5t(x) . (5.2) 

The implications of such variations on the Green functions Ga defined for operators of 
the form (1.2) with Dirichlet boundary conditions on 9A^ and other quantities, such as 
the extrinsic curvature Kij, whose definition also depends on the specification of 97W is 
considered in appendix C. The functional Schrodinger equation for \1/ has the general form 



(^'^+«°w)*(« = °' ^--""Wfe^' f^-^' 



where Tio is the corresponding Hamiltonian operator. At lowest order in the loop expan- 
sion, as shown subsequently, it is sufficient to let 

no^n = -y^h^ ^ + fi-'v{^) , (5.4) 

where V is a potential depending locally on ip and its derivatives in 9A1. Although \1/ is 
finite for arbitrary ip and smooth boundary surfaces 9A^ it is necessary to take account 
also of additional divergences arising from the second order functional derivative in Ti, 
in general 6'^/6ip{x.)Sip{x.') is singular as x' -^ x. In our calculations using dimensional 
regularisation is also sufficient to ensure a well defined Hamiltonian operator Tio rnay be 
found, satisfying (5.3), which differs from 7i in (5.4) by local counterterms which are just 
poles in e. These additional counterterms have essentially a similar form to the terms 
already present in H. As shown later the detailed structure of Hq is strongly constrained 
by the renormalisation group equations for \& expressing independence of the arbitrary 
scale //. 

26 



The semiclassical expansion to the wave functional ^((^) is derived by a loop ex- 
pansion of the functional integral (5.1) where the field (f) is also expanded about a back- 
ground classical solution ipc, with boundary conditions ^c\dM = <^- If the classical action 
5'c((^) = S{lPcj9), since the boundary term 5" in (2.1) vanishes for Lpc, then to one loop 
order the semiclassical approximation to \1/ is, from (2.10) with \l/((^) — expF^ip; (fc), given 
by 

^s.c.('^) = (det Areg)"^ e-*^"'^^(^) , (5.5) 

where v'c is the classical solution, with boundary value (^, and det A^g is the dimensional 
regularised functional determinant for Dirichlet boundary conditions defined as in (C.18). 
In order to discuss the Schrodinger equation (5.3) it is necessary to determine the effect 
on the classical solution (fc of independent variations in the boundary value ip and also of 
deformations of 9A^ as in (5.2). In both cases the variation d(fc satisfies A5(/7c = 0. For 
changes in (^ then the boundary condition is simply SipddM = ^^ and hence 

^'■Pci\X) , /\^/ I /r- cN 

with Gij{x,x') the Dirichlet Green function for A. For changes in the boundary dAi, as 
represented by (5.2), then we require 



and thus 



(v?c + 5v?c)L^,(^)+5,(^)=^(x) ^ 5v?c(^)L^,(^)=5t(x)a,V'e(x) , (5.7) 

^^ = G,j{x,x')t'r,\^,^^^^^dn^cj{^) . (5.8) 

From (5.7) and (5.8) it is easy to see that 

[-^-dn^ciyrj^cix) =0 , (5.9) 

which is important in later calculations. 

The verification that the semiclassical form for the wave function in (5.5) is an approx- 
imate solution of the Schrodinger equation, which becomes exact for quadratic Hamiltoni- 
ans, is well known in conventional quantum mechanics [20]. Here we extend the analysis to 
the next leading, or two loop, approximation which is straightforward but also pay special 
attention to the divergences which arise in a field theory context. To start it is necessary 
to consider how the classical action Sc depends on (p and also its dependence on changes 
in 9A1 as given by (5.2). By integration by parts and using the classical equation for t/^c, 
derived by requiring Sc to be stationary at (/> = (^, 

Sc{(fi) = dv {V{ipc) - lV'i{lfc)(fci) -\ I dSlfidn^ci ■ 

JM JdM 

With this result then under arbitrary variations in ip and also (/Jc, subject to /S.5(fic = 0, 
we find 



SSc{(p) = -^ / dS {6ipci + Sifi)dn(pc 

JdM 



'dM 
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Hence, taking Sipc = ^'P^ 

r 

——s,{^) = -a^v'c^(x) , (5.10) 

d(/7i(x) 

while, using (5.7) and also (C.8a,b) with (C.5) and the equation of motion for ipc on the 
boundary, we derive a local version of the Hamilton-Jacobi equation 

Sc{'p) = -ldn^ci{^)dnipci{^) + V{(fi{x.)) + ^Y'^dp(fii{±)dq(fii{x.) . (5.11) 



6t{x.) 



The essential results (5.10,11) ensure that at lowest order in (5.3) it is sufficient to 
take for V in (5.4) 

V(<^) = V{P) + ir'dp'fid.'f, , (5.12) 



S _ ,\ _i,,-iQf,%\ 1, S „ , ^, _i,,~eQf,^\ s 



since then 

(n^+7^)e-^^"^^^(^) = in^5e(^)e-*^'^^^(^), ^^5,(^) = -/C,,(x,x), (5.13) 

from (5.6) and (5.10) where /C is defined as in (C.22), 

/C,,(x,x') = 5nG',,(x,x')^'n|^^^(^)^,,^^(^,) . (5.14) 

From (5.9) and (C.21) we may find 

(^-9^</'c^^)lndetAreg = ^5c('^) + Y^tr(i34) , (5.15) 



where the pole term in e is necessary for the subtraction of divergences arising from b'^ jb(p^ 
in this case. Explicitly from (C.17), considering only terms which depend on V and are at 
most linear in K, and using the equation of motion on dM. Vn^c — Kdn^c+y^ = V'i'P), 

Buj = Iviivi:^ - \v;^vi; - ivV/' + iKdr^Vl'^ + . . . , 

tr {B^) = ^v:;{0)v:;i0) - ^v::;^{mji0) - \9^Jkk^n^c^^^^cJ (s.ie) 

which is a local function of ip and dn^Pc on dM.. In consequence from (5.13) and (5.15) 



\ ot / Ibn^e 



(5-17) 

60^ -^ -- -res 2 ^ -^ -- --reg ^ ^^^ -- ^^eg ^ 






Xi = \^^[ Y^ In det A^g — h -zA In det A^pp- ^r^ In det A 
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To first order in h Xx can be neglected and so ^s.c. satisfies the Schrodinger equation (5.3) 
in the leading semiclassical approximation if 



J 



7i 



(5.18) 
since, with these expressions, it is easy to see that 

1 A 

^2 = ^ 16^ (SijkkdnV^cz - ^K V-lj{0)) — In det A^eg , (5.19) 

At the next order in the semiclassical expansion of the wave functional we may write 

^{0) = ^,.e.(^) e^^^e^+o(a^) rS = r(2) - ^(2) , (5.20) 

where T^^) = r« + T^ + T^ + T^ as given in (2.13), (2.15), (2.17) and (2.18) while S^^\ipc) 
is the overall local two loop counterterm which subtracts the divergent poles in e. The 
calculations of section 2 determine the specific boundary terms. From (2.7,9) and using 
the equations of motion for (fc this may be expressed in the convenient form 

S^^\ip,)= f dvV^^\lJ--^'^,)+ f dS {-IJ--^'0^^^{^)dr,iPcz+ IJ^-'P^^H0)K) + . . . , 
Jm JdM 

(5.21) 

where Vo((/>o) = Ai~'^(^(</') + V'c.t. (</*)) ^'^d '^o ^^ given in terms of Z^'^^ and (p)^l^ by (2.25). 
In order to evaluate Xi, X2 in (5.16,18) we need to obtain an analogous formula to (C.19) 
expressing the dependence of det A^eg on (p. Using (5.6) and also the explicit form of the 
counterterms given by (1.17) in this case 



5 



In det Av 



(5.22) 
where Gmn\ is defined by (2.14) and /C is as in (5.14). Furthermore since 6G = —G6AG 

A r 

Gj^j^{xi,X2) = - dvGj^j{xi,x)Vjke{Lpc)GkJ2{x,X2)G£i{x,x)'dn\^^^^^^ , 



S(pi{x 
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we may also obtain in turn from (5.22) (as usual with dimensional regularisation we set 
5^~^(0) = so that the K term in (5.22) does not contribute to the second functional 
derivative here) 

in det Ar 



CIV (Jmn I y-^) Qmnij ^ jk\'^i •^) (JnOji^ki [p^^ '^J | £=2:fx) 

dvdv' Grnn\{x)Vll^^j{lfc)Gjj'{x,x')V-l'j^,t,,{lf'c)Gk'i{x\x)^ndn^ 

dvdv' dnGij (£, x)V-'^^{ipc) 

X G,r,^>{x,x')Gnn'{x,x')V^,^,^,{ip'c)Gj>i{x\x) 9n|^^^(^) 
^T67r% / '^'^^'^'^^^^^'^^^i'^r^{^''^y!!!nk{^c)Gk^{x,x)'^n\^^^^-^^ 

dS'dv Vl',.,j,{0')d'r,Gj'k{x', x)Vl.g^{ipc)Gmi{x, x)\dnGu{x, x) 



^g 2 9kkijlCji{k,±) . (5.23) 



In order to verify that (5.20) is a valid approximation to two loop order it is sufficient 

to show 

5 ^ J_ 

since this entails 
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(- - a„</'c^^) r(2) = -Xi -X2- A'loc , (5.24) 



{n- + n + n\}~^)i>{ip) = -{n''x,,, + o{K'))i>{ip) . (5.25) 

After using the equation of motion for ipc to eliminate Vn^pc we assume that Aioc(<^, dn'fic) 
is local function, which for the renormalisable theory considered here is of dimension 4. 
This may then be cancelled by an appropriate choice for 

v!i^ = n'x^oci^^p) , p = fi'h— , (5.26) 

dip 

by virtue of (5.10). We therefore endeavour to demonstrate (5.24) by considering each 
contribution to T^'^^ in turn, with ip — * ipc, and using (5.9) with (C.8a,b), (C.ll). The 
surface terms arising according to (C.8a) from the differentiation of the volume integrals 
expressing ra,r5,rc vanish since G satisfies Dirichlet boundary conditions. Beginning 
with the expression for Fa in (2.13) it is simple to obtain 

dn^Pcii^) 



5^(x) SP^{^)' .g27) 

dv gijk£Gij\{x) Gkm{x,x) dndnGmeix,x)\^^^^^^ . 
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For Ffe given by (2.15) the corresponding result is 



— / / dvdv'Vl!j'f,{ipc)Gii'{x,x')Gjj'{x,x') 

1 1 



(5.28) 



2 167r2£ / ^^^/jfc('^c)^»"£('/'c)G'fc^(x,x)a„anG'^£(£,x)|,^^^^^ . 

For the one particle reducible amplitude given by (2.17) we may obtain in this case 
5 ^ ... d 



(j^-^-^-^" 



.5t(x) 5(^i(x) 

^ ~R ^^^^'<5jfc|(x)y/Jfc(v7c)Gi^(2^,£)Ui9nGft'(£,xO^//j,fc,((/7c)<5j'fc,|(xO|^^ 



,s 



+ ^ j j dvdv'G,k\ix)V^'-kMGu'{x,x') 

X "^*'j''fc'(¥'c)Gj'^(a^', x) dndnGek'ix, x')\.^^^^^ . (5.29) 

For the remaining term in (2.18), using now (C.23) and (C.24) as well as (C.ll) and also 
Gm\ = 2ij-^Vkei0)/{16n^e) on dM. 



(^i)'^^'''^^^"^^)'^' 



1 1 



-gijkkdnVcjiy^) j dv' dr,G^i>{x,x')\^^^^^^ Vl!'j,f.,{ifi'c)Gj>k'\{x') 



4 167r% 

+ J dS"dv'V^';,i0")dnGjA£'\x')V^'j'k'i^c) 

X Gj'i[X ,X) OnOnGik'{x,X )\^^^(^^\^//^^/^,/\ 

+ 1 ds'ds" vi;'^{^")]cur , x)/CfcHx, ^)vi,],^,m 



Combining the results for the individual contributions we therefore find 

(^ - 9,V'c^^) r(2) = -Xi -X2-Xs- ^loc,l , (5.31) 
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where, neglecting K^ terms, 

Li ^ r 

- \Vl!j',{^)VZ{^)Vu{^) (5-32) 

+ l9ikkmgj£emY'^dp0^dq0j + lV-l!k{0)V--{0)V-u{0) 



This result is the critical part of verifying (5.24) but in addition we need to consider the 
contribution arising from the overall two loop counterterms as in (5.21) 

+ ^-i-y/(^)^i? - f,-^{p^^'{0)d^^,,K + V'p^'\0)) + . . . . 
Hence (5.31) and (5.33) together ensure the validity of the crucial equations (5.24) or 

(2) 

(5.25) with A'loc = '^loc,! + ^^100,2 which may then be used to determine Hq , as described 
in (5.26), so that ^((^) satisfies (5.3) at two loop order. 

Together with (5.17) the result of the calculations here are compatible with the 
general form 

Ho = - ^n'^ —^ -UK pQ^{ipo)j^— + Vo{^o) , , ^^, 

0^0 Oifoi (5.34) 

Vo((^o) = Vo((^o) + lY'^dpipoidqipoi - V'^po{0o) + • • • , 

for poi'po) = IJ'~^ pi.^) and where the neglected terms represent contributions to Vq involving 
scalars formed from the metric of dimension 2 or 4. The essentially complete expression 
for 7^0 to all orders, in terms of bare quantities, shown in (5.34) is a direct result of 
the requirement of consistency of the Schrodinger equation (5.3) with the renormalisation 
group equation (2.27) for the wave functional \1/ which is here written as 

Vi>{^) = i(C(/3^) +C(/3\/3^))^(^) . (5.35) 

Applying the renormalisation group operator "D, as given by (2.26), to (5.3) then since 
from (2.29a,b) VipQ = 0, Vpo = —p~^(3P we may use 

[V, ^] = , [V, n^] = hp-'K—-(3P{^) —- + fi-'V^(3P{^) + . . . , (5.36) 

and also, with C defined by (2.24), 

6 



[nl,^cm] = -^-^vV(^) + . 



to show how the form of Hq is dictated by imposing (5.35). The simple structure of Hq 
is at least in part a consequence of the restrictions on counterterms obtained by using 
dimensional regularisation with minimal subtraction. 
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6 Conclusion 

In this paper we have been primarily concerned with calculational techniques for 
quantum field theories with a boundary in the simplest case of a scalar field theory. The 
necessary labour is significantly increased beyond that for field theories without boundary. 
Even if the manifold on which the quantum field theory is defined is flat, as is usual 
for applications in statistical physics, the machinery associated with covariant treatments 
on curved space are a necessary ingredient in our method of calculating divergences for 
general curved boundaries. In spite of this complexity there are nevertheless potential 
formal developments associated with the functional Schrodinger equation discussed in the 
previous section. This may offer a non perturbative approach to discussing quantum field 
theories, such as in variational methods applied to the effective potential, but it is clearly 
important to ensure that the functional Hamiltonian operator is well understood in a 
perturbative context first. A further area of possible relevance of these calculations is to 
critical phenomena in real systems with boundary, as was the motivation for the work of 
Diehl and Dietrich [5] whose results for plane boundaries are reproduced here. Even for 
simple systems there is a rich variety of phase transitions associated with boundaries [4] 
whose universal properties may be calculated using renormalisable quantum field theories. 
For the 0{n) symmetric model discussed in sections 2 and 3 there is a fixed point for 
g = g^, m? = and on the boundary in the Dirichlet case (p = while in the Neumann 
case it is necessary to require (3c = or c = 0{K) (for a plane boundary c = 0). More 
generally if there is a cubic coupling gijk in Q{(p) in the treatment of section 3 then from 
(3.18) to one loop 

Ajfc = ~2^9ijk + if^„2 \^9iem9jmn9kn£ + ^9(.m{ij9k)m(. ~ 29Um{i9jk)m) ' (6-1) 



1 



SO that there is a possible fixed point for g^ = 0{e^), as discussed by Diehl and Ciach [22] 
for a single component field. 

In connection with the discussion in section 5 it is interesting to note that recently 
Liischer et al [23] have considered the Schrodinger wave functional for gauge theories with 
a view to applying finite size techniques to lattice gauge calculations. Although in other 
respects more complicated, due to requirements of gauge fixing and introducing ghost fields, 
the gauge theory case is simpler to the extent that there are no gauge invariant operators 
with dimensions < 4 which can occur as boundary counterterms. Lattice calculations of 
wave functionals are of course well defined non perturbatively but perturbative results, 
such as obtained here, provide useful constraints on the approach to the continuum limit. 



One of us (DMM) would like to thank the Royal Commisssion for the Exhibition of 
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Appendix A 

The short distance behaviour of the Green function Ga near the boundary may be 
found from an asymptotic expansion of the heat kernel Qa since 

/•OO 

Ga{x,x') = / dTQA{x,x';T) , (^-1) 

Jo 

where Qa satisfies 

(^— + A,)^A(x,a;';r) = 0, gA{x,x';0) = 5'^{x,x') , {A.2) 

with boundary conditions corresponding to (1.3) 
(l-P(x))^A(x,x';r)|^^^^^^=0, (P(x)n^(x)D^ + ^(x))^A(x,a:';r)L=,(^)=0. (A3) 

The extension of the DeWitt ansatz for Qa takes the form [13] 

^A(a^, x'; t) - ^ fe-"("'"')/2-0(a;, x'; r) + e-^("'"')/2"0(a:, x'; r)] , (A4) 

(47rT)2 V / 

where a{x, x') is the geodetic interval based on the geodesic path from x' to x while o"(a;, x') 
is the corresponding biscalar resulting from the geodesic which undergoes reflection on the 
boundary. The conventional DeWitt asymptotic expansion [14] is obtained by writing 
0(x, x'; t) ~ Aa (x, x') X]n=o '^n{x, x')r"', with Aa (x, x'), an{x, x') regular for x ^ x' and 
A2(a:,x) = ao{x,x) = 1. This expansion determines the singular behaviour of Ga{x,x') 
for x' -^ X, as in (1.5,6). The presence of the additional term O allows the boundary 
conditions (A. 3) to be satisfied by the full expression in (A. 4), while maintaining the usual 
form for O. An appropriate expansion of then enables the singular terms of Ga near 
the boundary to be obtained. 

To achieve this in the neighbourhood of the boundary, with coordinates so that the 
metric takes the form (1.1), both a and a are expanded in powers of e, according to (1.9), 
where to order e^ from I 

a = ^{y- y'f + \l^Ja'a^ - {y + y')\K,ja'a^ , (A.5a) 

^ = \{y + y'f + Wj^'&' - ^ ^^[ iK^j^'^' ■ (^-56) 

Using these results O may be obtained as a formal power series in e with 0„ = ©(e"^), 
taking also r = O(e^), by solving at each order sets of coupled differential equations in y 
and r with appropriate boundary conditions determined by 0„ in the analogous expansion 
of O. The solutions of these equations may be expressed in terms of the functions 

2r(?i+ l)/„+i(w, r) +ufn{u,T) = fn-i{u,T) + 2r5„o • 
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Using (A.5a,b) the results are more conveniently expressed in the form 

GAix, x'; r) ~ L_e--(-'- )/2- fe-"'/'^"0'(a;, x'; r) + e-"'/4"0'(a:, x'; r)] , (A.6) 

(47rr)2 ^ ^ 

where as before v = y — y' , u = y + y' . 

The first two orders in the expansion are straightforward to obtain in terms of the 
extrinsic curvature of the boundary Kij{-x) and also V'(x) and the projection operator 'P(x) 
appearing in the boundary condition (A. 3) 

n'o = i, vt\ = u\K,ja'an , n'Q = V-i, p_ = 2P-i, 

1 ^ .^ . (^-7) 

+ -U iK„aV^ V-I + /o(w, t){KV + 2^)/ 
r 

- 2P(7' A^ / + - Uu, t) {-y{l - P) A^ + y'Vb.Vy i . 
Using (A.l) and 

Jo (47rr)2 s! 

from the definition in (1.12), these results are sufficient to obtain (1.13). 

At order e^ the number of terms proliferate and the calculation becomes tedious. 
Here we consider only those terms which depend on X, which represents a mass^ term 
in the operator A, and also V', which has dimensions of mass, appearing in the boundary 
condition. For the conventional DeWitt expansion the relevant result is very simple 

0'2 = -rX°/, (A.8) 

where X^ = X\y=o. From the results of calculations carried out elsewhere we obtain 

0^2 = r(X° - 2VX^V)i - {y{\ - V)X^V + y'VX^il - V))fo{u, r) / 

+ iiyy' + 2r)/i(ti, r) - lOr^ fsiu, r)) {2K - ^K.^a'a^y I (A.9) 

+ 2r/i(w, r) {Ki; + 2^^) / + K.^a'a^i) I . 
The corresponding contributions to the Green function in the expansion (1.11) are 
G^=-\G^{v)X'i^ 

Gf = \G^{u) (X° - 2VX'^V)i - G'i,o(w) {y{l - V)X^V + y'VX'^d - V)) I 
+ G^i(ti)(|KV + 2V'') / + {Ayy'Go,i{u) + luGo,2{u))K,l, 1 
- {9,yy'G-i,i{u) + fCoaH + f wG'-i,2(tx) - GM)K^Ja'^^'^i . 
In (A.IO) 

G^{u) = G,{u) + ^ , GUu) = G,Au) + ^ ' 



(A.IO) 



which ensures that G|^, Gf'i have a smooth limit as £ ^ 0. 
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Appendix B 

In this appendix we describe the essential calculational details associated with de- 
termining the singular terms arising from the presence of a boundary, as £ ^ 0, of the two 
loop amplitudes considered in sections 2 and 3. For the graph in fig. la we need to consider 
the singular behaviour of the product of two coincident Green functions, i.e. G\{x) x G\{x), 
where G\{x) = G{x, x). From the results of appendix A the regular part of the coincident 
limit G\{x), as given by (1.5), has the general form near the boundary of the region on 
which G is defined 

G\{x)=G\{x) + ^{X{x)-\R{x)) 



2 



+ - (yr ^rf(x) - ^"'y+/i4(x)) + H{ 



X] 



e 



with Hj dyH regular at y = for (i ~ 4. y_|_ is a distribution or generalised function with 
support on y > and which has a simple pole at A = — p, p = 0, 1, 2, . . . [21], 



y+ 



^-1- ^ ^'^^^ 6^PHy) . {B.2) 



r-^ 



A + p p\ 
Furthermore it is important to note that 

- [vT-'-' - y'r'-') - 4 (- - ^) ^-^^''' (y) ' (^-3) 

with no simple pole in e, as may be verified by considering Fourier transforms. 
Writing for convenience 

as an expansion in e, and similarly for fd, Qd, hd, then from (B.2) and (B.3) the singu- 
lar contributions are, suppressing indices and assuming on the r.h.s. the products are 
symmetrised so that for instance cq x /o = /o x cq, 

G\xG\r^ |!^,' {-£go X eo^d"'{y)+eeo x f^ld"{y)-efo x fo^d'{y) 



+ 2{eQX go + e{ei x go - eo x gi))d'{y) 
-2(eo xho + foxgo 
+ e{ei X ho - eo X hi + fi X go - fo X gi 

+ ^ i-eoS'iy) + foS{y)) x H . (BA) 
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If we introduce the notation G\^ = G\y^o then from (1.6,7) and (B.l) we may write by 
analytic continuation in d from d < 2 

16^^ (S.5) 

In calculations of the divergent parts of amplitudes where G\ x G\ appear the non local 
singular terms present in (B.4) are cancelled by additional counterterms which involve 
contributions of the form 

-7^^{gox(^o+^°)5'(y)-(gox(/io + a^xO)+/ox(^o + ^°))5(?/)}. 

(S.6) 

Using the results from calculations of the asymptotic expansion of the heat kernel 
we obtain explicitly 

Co = "P- , ei = , 

/o = 4V' + |K + KV-, h = 4V' + |K + \KV- , 

^0 = -X^ + 2VX^V - 4V'^ - l^K , 

~g^= - (1 - V)X^V - VX^{1 -V)- 4V'^ - ^^K , 

/to = - dr,X^ + (1 - V)dr,X^V + VdnX^il -V)- 2P_XV - 2i)X^V- (S.7) 

+ SV'^ + ^V^K - 2PX°(1 - V)K - 2(1 - V)X^VK , 
hi= - \{dnX^ - 2VdnX'^V) + i(l - V)dnX'>V + \Vdr,X^{l - V) 

+ 2XV + S-^X" - f V'^i^ 

+ \X^K + VX^VK ~^[VX^{1-V) +VX^{1-V))K , 

where in ^o,i5 ^o,i we have given only those terms depending on X, t(j. For Dirichlet 
boundary conditions, V,ip = and from (B.5) 

G\^ = H^ = 0^ Q^{^G\f = dnH'' = Q, {B.8) 

so that the last line of (B.4) vanishes in this case. For Neumann boundary conditions 
P = 1 and 

dn{G\)^ = -V'G'I" - G|V , dnH^ = -ijH^ - i^V . (S.9) 

When integrating by parts in integrals over dv it is useful to note that 

dv = dSdy{l-yK + 0{y^)) , d'{y) =V ^{n^S{y)) + KS{y) , (S.IO) 
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with a corresponding formula for S"{y). 

For the amphtude corresponding to the graph in fig. lb it is necessary to determine 
the singular behaviour of 



J\S 



/ dydy' f{y)f{y')GA{x,x 

2^-' r V7 / ^^'"'^ 

- ,„ n I dy^f{yfil^GA{x,x) + IGa{x,x)1 + GA{x,x)l^)5dM{^,^') , 
IGyr-^e J ^/g 

where ^ = \^{l — yK + . . .) and the second term subtracts the subdivergence as usual. 
f{y) is an arbitrary smooth test function vanishing rapidly as y ^ oo and 5qm{'^t^) = 
^d-i^y^ — x')/a/7(x) denotes the surface 5 function on dM.. The singular terms as £ -^ 
arise for x ^ x' and it is convenient to use directly the representation of Ga provided by 
the heat kernel, as in (A.l), where writing 

is equivalent to a Taylor expansion of the Fourier transform. For the leading term in the 
expansion (1.11) it is necessary in (B.ll) to take Ga -^ Gi{v)+V-Gi{u) with Gi{v)\ = 0. 
Using (B.12) gives 

^''"'^ ~ '''yj^"' ^;(p^{r(d-|)IH-'^+%M + |r(d-|)|.|-"'+'v%^, + ... } . 

(B,13) 
dropping the R term. This is then easily shown, by considering specific test functions such 
as f{y) = e~py , to lead to a simple pole in £, 



JJdydy'f{y)f{y')Gi{v) 



r-w^ 



-2s 1 /.oo ,.00 (5-14) 

^(Te^i/o ^yf^y^'^'-l dyf'iy)"-mf'(o)ydM, 

where the singular contributions arise both from the usual short distance limit y -^ y', as 
for manifolds without boundary, and also y,y' — * 0. Similarly for Gi{u), although only 
the latter contribution arises in this case, 

dydy' f{y)f{y')Gi{uf r^ ^ —^ f{0)f'{0)6dM ■ {B.15) 

In order to discuss the Gi{v)'^Gi{u) contribution we write, using (B.12) again, 
Gi{vf-^6iv)5dM)Gi{u) 

li-lwT-^/ di ds ^_ ...,,,_,, e- -*- (5.16) 



647r^ ,^olr(d-2)y, io (s + t)^(<^-i) 



fj.- 



-e 



/oo 



e 
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where e plays the role of an infra red regulator so that for e > the r.h.s. of (B.16) has 
a rapid fall off for y,y' ^ oo and we may therefore take f{y) -^ a + by which allows all 
integrals to be explicitly evaluated. When e -^ the e dependence cancels in the terms 
containing poles in e, since they arise only from short distance effects, and we find 



^2 ^ %/ NX A^ / N ^ ^' 



dydy' f{y)f{y') (G^{v)' - |-^ d{v) ddM)Gi{u) ~ -^^^-^(2 - e) f{0)f\0)d9M • 

(S.17) 
For the final leading order contribution we obtain 

dydy'f{y)f{y')Gi{ufGi{v) ^ ^ ^^^ /(0)/'(0)5aM ■ (5-18) 

For pure Neumann, Dirichlet boundary conditions the surface contributions combine to 
give 

[ fdydy' fiy)f{y') ({G^iv)±G^{u)f T^S{v)^9MG^{u: 

^ .{±l-^e^J^e)f{0)f'{0)69M. 



(167r2£)2 

For the next order in the expansion in (1.11) we need to calculate similar results in- 
volving Gf and Gf as given in (1.13). Although the labour involved increases significantly 
similar methods may be applied, in particular when there is a subdivergence introducing 
the regulator e as in (B.16) so that we can take now f{y) -^ 1 and evaluate all integrals 
before letting e -^ 0. 

2 /^B 



JJdydy'f{y)f{y')Gi{ufG^ 



^ ^ j^, fiOniKV. + IKV + 2i;)59M , (S.20a) 

dydy'f{y)f{y')G^{u)G^{v)Gf 

IJdydy' f{y)f{y') [Gi{vf - £^ d{v) 5a^)Gf 



fi- 



-2e 



(167r%)2 
j jdydy'f{y)f{y')G,{ufGf 



f{Qf{\{l-%e)KV- + ^{l-\e)KV + Ai;)5dM , (5.20c) 



^ ^ _L_ /(o)2 IK SoM , {B.20d) 



e (167r2)2 ^ ^ ^ 4 
IJdydy' f{y)f{y')Gi{vfGf 
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^-'^ 1 



■{j dyyf{yf{KV^ + 2K'^V^Vj)- j dyyf'{yfKy9M 



6e (167r2)2 

(S.20e) 



jldydy'f{y)fiy') (Gi{v)Gf - J^yK 6{v) 6dM)Gi{u) 



^■ 



-2e 



(167r% 



2c-^2 



f(0fl{l-le)K59M. (S.20/) 



The terms on the r.h.s of (B.20e) correspond to the expansion of V'^S'^ around y,y' = to 
this order in accord with (1.8), d'^i^x^x') = (1 + yK + . . .)d{y — y')5a7V((x, x'). 

For the remaining one particle reducible graph in fig.lc it is necessary because of the 
form of the counterterms to consider the Neumann and Dirichlet cases separately. The 
basic Green function may be written as 

Ga(x, x') = (1 - P)G'^(?/, 2/0(1 -V)+ VC'iy, y')V + . . . , (S.21) 

where other contributions are irrelevant in the following discussion and the explicit depen- 
dence on X, x' is suppressed. In the Neumann case it is sufficient to investigate integrals 
over y,y' of the essential form, assuming Pd{,y) = f{y){Gd + yfd) where ed(x), /d(x) are 
given by (B.7) for P ^ 1 to the first order in an expansion in e and f{y) is a suitable test 
function as considered previously. 



rN 



jjdy dy' yl-" y J"" Pd{y) G'^iy, y') pa{y') 



-e 



^ ^ jdyyl-''pd{y)G''{y,0){pUO)-pmi;) 

+ fdy'y'^-''{p^iO)-p,{0)ij)G''iO,y')pd{y'] 



(S.22) 



+ ^ (p4'(0) - p4(0)^)G^(0, 0) (p4'(0) - P4(0)^) . 

As remarked in section 3 for the regular part of G^, using the appropriate boundary 
conditions, the counterterms in (B.22) may be integrated by parts to give 

^"^'"^ ^ IJdydy'Diy) G'^^^^^y^y') D{y') , D{y) = yl-'^p,[y) + ^5'{y)p,{y) , (S.23) 

where by virtue of (B.2) -D(y) is regular, regarded as a distribution on smooth test func- 
tions, as £ — * 0. In order to calculate the poles in e arising from this one particle reducible 
amplitude then it is necessary to only consider the leading singular contributions to G^ 
when y,y' ^ and also x' -^ x, as exhibited in (1.13) or (1.15). 

The most singular term in (1.13) corresponds to taking G^{y,y') -^ Gi{v) +Gi{u) 
and also letting -0 — > in (B.22). Analysing the integrals over y,y' in this case gives an 
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expansion in powers of a of the form 

- J -^(7(pd(0)p.'(0) + p,'(0)pd(0))(2a)i-i'^ (5.24) 

- 6{pMpm + p4'(0)pd(0))(2a)i-<^) + . . . , 
for 



27r^<^ 



The remaining terms not shown in (B.24) are less singular as a ^ and do not lead 
to poles in e. The numerical coefficients a, 7 and 5 may be easily evaluated in terms of 
gamma functions for general d and as e ^ 

a = 0(1), 7~l-£, 5 r^ ^-%l + 0{e^)) . (S.25) 

Hence the poles in e appear to cancel as expected by virtue of the counterterms present 
in (B.22). Nevertheless if in the expression obtained in (B.24) {2a)~^ is regarded as a 
distribution on dAA defined by analytic continuation in s from s < ^{d—1) then it is non 
singular when s ~ 2, as in the first term of (B.24), but for s ^ |((i — 1) there is a pole 
given by 



e 



((2a)-^(^-i)+i- - (2a)-i('^-i)+i^^) ~ ^ (^ - \)Sa-i ^M + 0(1) . 



This result is a direct extension of (B.2) and (B.3). Hence although I^ is a one particle 
reducible amplitude including counterterms for all sub- divergences it still results in a local 
divergence represented here by poles in e of the form 

(S.27) 
where we have inserted the explicit form of Pd{y), with f{y) an arbitrary smooth test 
function, and also used ei = from (B.7). For the terms proportional to t/j then in 
(B.22) from (1.13) we let G^{y,y') -^ 2i/jGi^o{u) or in the ifj dependent counterterms 
G^{y,0) -^ 2Gi{y), G^{0,y') -^ 2Gi{y'). In a similar fashion to the analysis leading to 
(B.24) and (B.27) 

4^ ~ - V' ^ (27 Pd(0)2(2a)i-i^ - 5{pa{0)p4{0) + pmpM) (2a) t-^) + . . . 

2e 

—^tljfiOfeoeoSdM, (S.28) 
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where we have used (B.26) again with 7 , ijl'^5 = 1 + 0(£^). The analogous results for the 
K dependent terms are obtained in (B.22) by taking 

C'iy, y') ^-^ {2cj{Go{v) + Go{u))u + 2GMuyy' + {d- 2)^,0^) 

^ (S.29) 

The various additional terms represented by X, which may easily be found from (1.15), 
may be discarded here as they do not result in any overall divergence since on using 

(2a)-^(^+i)+^^a.a, ~ \^1^J ^dM • (5-30) 

in conjunction with (B.26) ensures that the poles in e cancel. In this case the counterterms 
in (B.22) are also irrelevant in obtaining the leading singular term as a ^ and we obtain 



K 1 /n9/^n93j U ^ 



/" = -!-^^^K.pi(Qf(2a)i-i-' ~ ,jl , ieA'/(0)^e„e„feM , (B.31) 



since «; = 3 + 0(e) (individual contributions to k possess poles in e but these cancel when 
combined together). Combining this with (B.27) and (B.28) gives finally for the singular 
part of the integral in (B.22) /^ = /f + /f + /f . 

In the corresponding Dirichlet case we are concerned with integrals of the form 
/^ = lldydy'yl-'y't%a{y)G''iy,y')p,{y') 

jdy yl-'' pa{y) G^'iy, y')%'\y,^, P4(0) 



+ ldy'y'i-''p,{0)dyG''{y,y\^^p,{y'] 



(S.32) 



,, -2e 



+ ^ pM dyG'^iy, y') dy,\^^^,^^ P4(0) , 

where paiy) = f{y){^d + yfd) as before but now P ^ 0. The analysis of the singularities 
of I^ follows similarly to the previous discussion for I^ . Again for the regular part of 
G^ it is possible to write (B.32) in the analogous form to (B.23), assuming G^'^'^^^O, y') = 
G^'''^s(y, 0) = 0, showing that poles in e can only arise from the singular short distance 
part of G^ as can be obtained from (1.13) or (1.15). In this case it is necessary to note 
that dyG^ dy'\y=y'=o, wWch appears in the 1/e^ counterterm in (B.32), involves terms 
proportional to (2(t)~2<^^ (2a")~2(<^~i) when x' — ^ x, as shown in (C.25), which represent 
non integrable singularities on integration over dAi. Furthermore the resulting divergences 
are not regularised by analytic continuation in d. Elsewhere [17] we have discussed how 
to deal with this problem, which also occurs in the calculation of electrostatic energies on 
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Ai when the potential is given on dAi. It is sufficient to treat the singular short distance 
terms appearing in dyG^ dy'\y=yf=o as distributions but it is crucial if correct results are 
to be obtained to also include an additional term involving Sqm proportional to K which 
is shown in (C.25). 

For the contributions arising for G^{y,y') -^ Gi{v) — Gi{u) then instead of (B.24) 

+ -^(7(pi(0)Pd(0)+Pd(0)pi(0))(2a)i-i'^-5(pi(0)p4(0)+p4(0)pi(0))(2a)t-^ 
+ ..., (S.33) 

where the numerical coefficients here have the form 

a, 7~l-£, (3 ^ ^~^{l-\e) , 5 ^ ^-^{l-e) . (S.34) 

In this case the counterterms contribute to both the leading singular parts shown in (B.33). 
These are necessary to ensure that only poles in e with local residues may appear. Using 
equations (B.26) and (B.34), the result becomes, with the explicit form of paiy) again, 



la ~ - (i6^2g)2 (/(0)/'(0) (l-£)eogo + |/(0)'((l-£)(eo/o + /ogo)+£(go/i+/igo)))59M , 

(S.35) 
For the K dependent terms then using instead of (B.29) 

G^(2/, y') ^ ^ {2a{Go{v) - Go{u))u - 2G^{u)uyy') + X' [k - {d - l)^^^g^ 

yields, since no relevant contributions result from X' , in the Dirichlet case 

4'' = J^;5^(2^M0)^2a)i-i^-A(M0)p4(0)+p4(0)p,(0))(2a)i-^) 

9. ^ ^ (5-36) 

where the last term is is a consequence of the extra 5qm contribution shown in (C.25). 
The coefficients are «;, /x'^A ~ 1 — 2£ so that (B.36) gives 



-2e 

(167r"2F)^ 



Ih ~ ,,^ 2^2 ¥K /(O)' eo^o hM . (S.37) 



The total singular contribution resulting from (B.32) is then I^ = I^ + l/^ as given by 
(B.35) and (B.37). 
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Appendix C 

For the derivation of the functional Schrodinger equation in section 5 we here sum- 
marise the necessary results for variations induced by a deformation of the boundary 
surface dAi. For dAi parameterised by coordinates x^ and the embedding in terms of the 
coordinates for Ai specified by x^(x) we define 

e^,(x) = -^ , n^e^, = , n^ = 1 , (CI) 

with n'^(x) the unit inward normal on dAi. For a metric gr^;^ on At the induced metric on 
97W is determined by 

iC.2) 





lij 9fii^^ i^ j 1 


V.e^,- 


= d,e^, + r^^e^.e", = e%V% + n^K, 


V,n^ 


= -i^^,e^^ d, = ehd^ 



and then 

Y7.cM. — /^.cA*. _|_ FA* c'^.oP. — cA*. pfc _l_ ^A* fc 



(C.3) 



defines the extrinsic curvature Kij = Kji of 9A^ for F^ the usual Christoffel connection 
formed from g^^j. If Vi is defined with the connection F^- = F^^ acting on tensors on dAi 

then Vi^jk = so that F^ is the Christoffel connection formed from 7j/j. From (C.3) it 
is straightforward to derive the Gauss-Codazzi equations relating the Riemann curvature 
tensor R^uap for x G 9A^, with zero and one component along the normal n^, to the 
intrinsic Riemann curvature Rijue. of 9A^ associated with the covariant derivative V^ and 
also the extrinsic curvature Kij. 

For a variation 5x^ = —n^5t, as in (5.1), we let 



5tK,je^^ ~n^di6t 



(CA) 



from (C.3) and hence, regarding n^ as defined by (C.l), 

S'n^ = 6n^ - n,T';^n''6t = e^d.St . (Cb) 

With the induced metric given by (C.2) we have 

6% = 2K,j6t, {C.6) 

and using Kij = n^Vie^j from (C.3) 

For integrals over a local scalar function / on A^ or, restricted to the boundary, on 9A1 

S f dvf= [ dSStfldM , (C.8a) 

Jm JdM 

S [ dSf\9M= [ dSdt{-dnf + Kf)\ (C.Sb) 

JdM JdM 
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neglecting in (C.8b) any dependence of / on n^, 7^^, Kij . . . where it would be necessary 
to use results such as (C.4,5,6,7) above. 

The deformation of the boundary induces a variation in the Green function Ga and 
the heat kernel Qa since boundary conditions are an essential part of their definitions in 
(1.4) and (A. 2). In the Dirichlet case, with a choice of gauge so that A^ = as in (1.10), 

{SGAix, X\ t) - 6t{x.) dnGA{x, X \ t)) |^^^(^)= , (C9) 

and hence for arbitrary x^x' ^ M. 

6Qa{x,x']t) = / dr' / dS 6t{k) gA{x,x;T - r') dndnGA{x,x';T')\ ,.. . (CIO) 
^0 JdM ^ ' 

Analogous results hold for Ga as a consequence of (A.l) and therefore 

-t-7TtGa{x,x') = G'A(a^,^)5n5nGA(£, a^')L=^(x) • (Cll) 

Since by definition Qa is the kernel for the operator e~'^^ it follows from (C.IO) that 

5TV(e-^^)=r / dS5t{yi)ix[^^gA{x,x'■T)"^'r.l^^,,^^ . {C.12) 

JdM ^ ' 

In general the asymptotic form of the heat kernel is of the form 

(47rr)^"Tr(e-^^) ^ \ dv VS2n(x)r-+ / dS Vs^(x)r^" , (CIS) 

Jm ^=0 JQM n=l 

and we may also write 

T{^'KT)^''dr.QA{x,x'-T)"d'r\^^^,^^^^^ ~ ^i3.(x)r^- , (^.14) 

n=0 

defining i3n(x) on dM.. As a consequence of (C.12) 

Bn + ^^ [ dSB^ =triB^) , (C.15) 

^t JdM 

with Bn = for n odd. This result may be verified using our previous results [13] for 



n = 0,1,2, 



4 = 1, Bi = -\V^K , B2 = \{R + K^- K^'K.j) - X\dM , (C.16) 

with R the scalar curvature on dAi . For the discussion in section 5 we need the contribution 
to B4 given by 

i§4 = {^x^ - Iv^x - \b^x + iKdnX) |g_^+ . . . , (ai7) 
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where on dAi V^f X = n^n^'^ yV^X. This has been determined by direct calculation and 
also checked by using (C.15). 

For the functional determinant defined by (1.15) the finite regularised form in four 
dimensions is given by 

lndetAreK = lndetA + ^^| / dvB^+ I dS bA , (C.18) 

167r^£ yjM JdM J 

where S4, S4 are given by (1.16). For general boundary conditions, with A given by (1.2), 
it is easy to see that, under smooth variations in X, 

51ndetAreg= / dvtr{DASX), {C.19) 

Jm 

where Da{x)j as given by (B.l), (B.7) and (B.12), has the form 

Da{x) = Ga{x,x) + -^{X{x) - ^Rix)) 

ibw^e ^^20) 

+ -[^,(;^f^i^'i^)'P-i^)^n{x))-{4i;{ic) + |K(x))5.(x)) , 

with Sn{x) the surface S function defined by dv{x)dn{x) = (i5(x), in the standard coor- 
dinates of (1.1) Sn{x) = S{y), V^(n^(x)5n(x))= S'{y) — K{x)5{y). This represents the 
regularised form of Ga{x,x) in dimensional regularisation with both short distance and 
boundary divergences removed by subtraction of poles in e. In the Dirichlet case, when 
V- -^ 1, we may also find for variations of the boundary, from (C.12) and (C.15), 

^ Indet Areg = -tr (/Ca(x,x)) + :^ tr (i34(x)) , (a21) 



5t(x) ^^^ ^ ^^ ' ^^ ' 167r2£ 

where formally 

/Ca(x,x') = 5nGA(a:,a:')V„L^,(^),,,^,(^,) ■ (^22) 

For applications in section 5 it is necessary to also consider variations of the Green 
function when arguments lie on the boundary. By careful analysis, assuming /(x) and 
f{x) are smooth test functions on dM. and Ai, 

^ jdS'dv" /(x')5nG'A(x', X")\^,^^^^,^f{x") 

= dS'dv" f{k')}CA{y^',i)dnGA{x,x")\^^^^^^f{x") ^ ■ '' 

+ /(x)/(a:(x)). 
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Also, in a similar fashion 
6 



dS'dS" /(x')/Ca(x', x")/(x' 



5t(x) 

"dS'dS" f{Sc')ICA{Sc', x) /Ca(x, x")/(x") ^'-^■^^^ 

-r^(x)a./(x)a,/(x) - /(x)X(a;(x))/(x) , 

which is analogous to a result quoted by Symanzik*. 

In general the coefficients Bn are related to the short distance expansion of /Ca(x, x') 
when x' ^ x. We obtain 



JCA-^{2a)-'^'l+4^K^^ 



7r2 



id 



2{d-l] 



dM 



(C.25) 



The term proportional to 5qm results from a careful treatment of the singular behaviour 
of /C at coincident points [17]. The result given is compatible with treating {2a)~^ as 
a distribution [21] or in integrating the kernel /C(x, x') with a smooth function /(x') 
excluding a disc of radius c centred on x and subtracting a term oc c~^ to ensure a finite 
result as c ^ 0. 



See eq. (6.6c) in ref. [10] 
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